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ABSTRACT 


A  multi-stream  diffusion  model  is  proposed  for  the  calculation  of 
the  annealing  behavior  of  boron  that  is  ion  implanted  into  silicon  at 
room  temperature  and  subsequently  annealed  .  This  model  is  capable  of 
predicting  both  the  redistribution  and  the  electrical  activation  of 
boron  during  the  anneal,  as  a  realistic  model  should.  The  calculated 
results  compare  very  well  with  extensive  experimental  data  reported  in 
the  literature.  The  comparison  includes  samples  that  are  implanted  at 

room  temperature  with  boron  in  the  dose  range  from  1C)44  to  1046  ions/ 

2  . 
cm  and  subsequently  annealed  in  the  temperature  range  from  800  C  to 

1000°C.  This  range  of  dose  and  annealing  conditions  includes  both  the 
typical  applications  of  ion  implantation  as  it  is  applied  in  the  fab¬ 
rication  of  devices  and  the  unconventional  cases  of  high  dose  implants 
and  low  temperature  annealing. 

The  annealing  model  is  in  essence  a  diffusion  model  extended  to 
include  the  following  situations: 

(1)  precipitation  of  boron,  when  the  boron  concentration 
level  exceeds  the  solid  solubility  limit,  in  high 
dose  implantation  cases 

(2)  trapping  of  boron  by  the  strain  fields  associated 
with  dislocation  dipoles  that  form  during  annealing 
at  low  temperature 

In  both  cases,  the  presence  of  this  immobile  and  electrically  inactive 
boron  alters  significantly  the  redistribution  and  electrical  activation 
behavior  of  boron  during  the  anneal. 

The  presence  of  boron  precipitates  makes  it  necessary  to  include 
four  species  in  the  annealing  model.  These  species  are  substitutional 
(electrically  active)  boron,  boron-vacancy  pairs  (electrically  inac¬ 
tive),  positively  charged  vacancies,  and  the  aforementioned  immobile 
boron.  The  electrically  active  boron  is  assumed  to  diffuse  substitu- 
tionally  by  means  of  random  encounters  with  neutral  vacancies.  The 
boron-vacancy  pair  is  assumed  to  diffuse  much  more  rapidly.  The  prin¬ 
cipal  interaction  among  the  species  is  the  reaction  of  active  boron 
with  positive  vacancies  to  form  BV-pairs. 


The  parameters  in  the  model  art  a  set  of  dif fusivities  and  lifetimes 
for  various  species  and  interaction.  These  parameters  are  estimated  from 
the  examination  of  each  diffusion  and  each  interaction.  The  values  of 
the  most  important  parameters  obey  simple  activation  energy  relations.  In 
other  cases,  the  modelling  was  either  not  carried  to  the  extent  necessary 
as  to  make  the  temperature  dependence  apparent,  or  else  the  parameters 
represent  the  composite  effects  of  very  complex  interactions.  In  all 
cases  ,  the  mathematical  model  and  the  parameter  set  are  universal  in  the 
sense  that  the  same  equations  and  parameters  can  be  used  for  the  predic¬ 
tion  of  ordinary  diffusion,  proton  enhanced  diffusion,  and  the  annealing 
behavior  of  ion  implanted  boron  at  room  temperature  under  conditions  that 
include  precipitation  effects  and  low  annealing  temperature  anomalies. 
Furthermore,  in  most  cases  of  practical  importance,  high  dose  and  low 
annealing  temperature  conditions  are  of  limited  interest  and  ,  as  a  result, 
the  three  stream  diffusion  model  with  an  appropriate  subset  of  parameters 
is  sufficient.  This  reduced  model  is  still  very  important  and  attractive 
because,  once  the  parameters  associated  with  high  dose  and/or  low  temper¬ 
ature  anomalies  are  deleted,  the  remaining  parameters  have  simple  activa¬ 
tion  energies  and  are  determined  by  the  specification  of  a  single  varia¬ 
ble:  the  temperature. 
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The  temperature  dependence  of  parameters 


Furthermore,  it  is  possible  that  ordinary  diffusion,  proton  enhanced 
diffusion,  and  the  annealing  of  boron  are  indeed  the  same  problem  with 
different  excitation  and  initial  condition.  With  this  preliminary  dis¬ 
cussion,  we  turn  the  attention  back  to  the  consequences  of  a  perturba¬ 
tion  of  the  boron-silicon  system  at  equilibrium  with  proton  bombardment. 

Fladda  et  al  [111  have  shown  that  substitutional  boron  in 
a  diffused  sample  can  be  rapidly  removed  from  their  lattice  sites  by 
proton  bombardment  at  room  temperature.  This  same  result  can  be  achieved 
if  the  diffused  sample  is  replaced  by  an  ion  implanted  boron  doped  sample 
that  was  fully  annealed.  However,  if  an  unannealed  ion-implanted  sample 

is  used,  no  change  in  the  fraction  of  substitutional  boron  can  be  detected 

17  ,  2 

even  for  proton  doses  as  high  as  10  protons/cm  .  Similar  results  can 
be  obtained  at  higher  temperatures.  Figure  3  shows  profiles  from  a  boron 
diffused  sample  irradiated  by  10  KcV  protons  at  550 °C  for  two  hour4--  ,  ns 
determined  by  SIMS  and  electrical  resistivity  measurement  [12].  The  im¬ 
purity  distribution  is  characterized  by  the  presence  of  a  dip  in  the 
electrical  carrier  profile  near  the  proton  range  and  by  some  diffusion 
of  the  tail  of  the  profile.  This  dip  is  not  observed  in  tin'  total  boron 
concentration.  The  enhanced  diffusion  of  the  tail  is  more  pronounced  at 
higher  temperatures  such  as  850°0  02’  900  °c .  If  the  proton  dose  is  in¬ 
creased,  then  a  peak  starts  to  emerge  from  the  bottom  of  the  dip.  This 
peak  is  attributed  [13 1  to  pileup  of  immobile  boron  (complexes,  precip¬ 
itates,  dislocation  trapping,  etc.).  To  elucidate  this  issue,  we  analyze 

the  annealing  of  damage  produced  by  the  room  t empora turo  irrad iation  of 
1 5  2 

a  5  v  10  protons/cm  dose  on  an  annealed  sample.  The  transmiss ion 
electron  micrograph  shows  the  formation  of  linonted  defect',  that  we 
Identify  as  dislocation  dipoles  [191. 

c.  Ih’o  filing  Measurements 

Figure  1  shows  a  collection  of  isochrona 1 Iv  annealed  boron 

profiles  (35  min)  at  Non  V  to  10M0°(:  and  the  correspond  i  ng  electrics! 

*  11  10 

carrier  profiles  for  implantation  doses  of  in  to  in  ions. -cm  ,  after 


The  calculation  based  on  the  integration  of  the  SIMS  prof i le«  shows  that 
the  actual  implanted  do-os  are  3<>  ‘  f  ~  1(1',  and  30"  higher  than  1  he 
nominal  dimes  of  in'  in'--’,  and  l<'lo  ions  'em, ~ . 

1  1 


system ,  tho  equilibrium  condition  can  be  altered,  by  an  excess  vacancy 
concentration  produced  by  proton  irradiation.  Since  the  production  of 
vacancies  that  result  from  collision  is  a  complicated  problem,  it  is 
pertinent  to  elaborate  upon  the  conditions  under  which  excess  vacancies 
are  produced.  Obviously,  we  want  to  avoid  producing  severe  damage  to 
the  crystal,  since  in  this  case  the  recovery  of  equilibrium  will  be  a 
much  more  complicated  matter. 

Energetic  light  ions  such  as  protons  or  boron  produce 
diffuse  disorder  clusters.  These  clusters  are  isolated  because,  even 
at  room,  temperature,  there  is  a  certain  amount  of  annealing  that  pre¬ 
vent  the  overlapping  of  individual  disorder  clusters.  In  contrast,  if 
flu*  implantation  is  Der formed  at  liquid  nitrogen  temperature,  then  the 
disorder  clusters  will  overlap  and  an  amorphous  layer  will  form.  This 
i-  a  direct  consequence  of  the  negligible  annealing  at  low  temperature' . 
An  amorphous  layer  is  also  observed  if  the  implanted  ion  is  heavy.  In 
this  case,  even  if  the  implantation  is  carried  out  at  room  temperature 
mu.  ;  .<■  implanted  dose  is  moderate,  the  si/e  of  the  disorder  clusters 
are  so  extended  that  they  will  overlap. 

In  a  Proton  Enhanced  Diffusion  (PED),  generally  the  sub¬ 
strate  is  held  at  an  elevated  temperature  and  the  diffusion  occurs  under 
a  steady  proton  bombardment.  In  view  of  this  continuous  excitation  and 
tin*  criterion  of  thermal  equilibrium,  we  can  classify  the  PED  as  a  non- 
equilibrium-steady  state  problem.  We  shall  see  that  the  annealing  of 
ion-implanted  boron  bears  some  resemblance  and  some  differences  with 
this  case.  As  in  the  proton  case,  vacancies  produced  by  the  light  boron 
damage  give  rise  to  an  initial  nonequilibrium  condition;  and,  in  con¬ 
tra  -t  with  the  PED  case,  the  perturbation  is  not  at  steady  state.  Hence, 
annealing  of  ion-implanted  boron  is  a  transient  problem  with  a  nonequ i - 
Librium  initial  condition.  In  this  analysis,  the  simplest  case  is  tho 
ordinary  diffusion  occurring  under  thermal  equilibrium,  then  in  Increas¬ 
ing  order  of  difficulty:  PED,  a  steady  state  nonequilibrium  problem , 
and  finally  the  annealing  of  ion-implanted  Boron,  a  transient  problem. 

This  analysis  provides  a  natural  way  of  conceiving  how 
the  enhanced  diffusion  at  tho  outset  of  the  anneal  can  turn  into  ordi¬ 
nary  diffusion  at  a  Later  time.  In  thi--  interpretation,  it  is  simply 
a  transient  problem  reducing  t>>  the  thermal  equilibrium  condition. 


lo 


T 


T 


T 


"71 


Table  1 


A  COMPARISON  OF  DIFFUSION  COEFFICIENTS  AT  VARIOUS  TEMPERATURES 
OBTAINED  UNDER  DIFFERENT  CONDITIONS,  AFTER  HOFKER  ET  AL  [Hi 


D  (cm^  s 

Anneal  - 

N(|,0):  Not  Annealed 

ing 

Distribution 

N  (  d  ,  0  )  : 

1 o ni  0  • 

(  °r'\ 

Annea led 

from  Kurt/. 

\  ^ ) 

,  .  dose  : 

dose  : 

Distribution 

in14  . 

10  10ns  cm 

10l°  ions  cm“^ 

-15 

-16 

800 

— 

— 

1-5  .10 

5-2  .10 

(extra pol  .  ) 

-1-1 

-14 

-15 

-15 

900 

1-7  .  10 

1*4  .10 

1-4  .  10 

1-3  .  10 

(extra pol  . ) 

,  -11 

-14 

-14 

-14 

1000 

4.8  .  10 

5.2  .  10 

1-8  .10 

1.7  .10 

-13 

-13 

-13 

1100 

1.9.  10 

1-5  .  10 

1.7  .10 

b  .  Proton  Irradiation 

In  the  preceding  review,  we  have  examined  situations  in 
which  the  boron-silicon  system  is  near  thermal  equilibrium  conditions. 

In  other  words,  equilibrium  has  been  achieved  with  respect  to  tempera¬ 
ture  and  all  possible  chemical  reactions  in  the  system.  We  purposely 
exclude  the  attainment  of  equilibrium  with  respect  to  the  distribution 
of  matter  within  the  system  in  this  definition.  In  what  follows,  wo 
turn  the  attention  to  the  perturbation  of  the  equilibrium  condition. 
First,  we  discuss  the  use  of  protons  to  create  an  excess  vacancy  con¬ 
centration,  and  then  we  compare  the  differences  between  a  proton  en¬ 
hanced  diffusion  and  the  annealing  problem  with  respect  to  thermal 
equil ibrium . 

The  purpose  of  perturbing  a  system  at  equilibrium  with  an 
excitation  is  to  observe  the  subsequent  response.  Generally,  .1  system 
under  perturbed  condition  will  expose  the  driving  forces  that  will  ena¬ 
ble  the  system  to  return  to  equ  i  1  ibr  ium ,  hence  elucidating  physical  •••oeh- 
anisms  hidden  under  thermal  equilibrium  conditions.  In  the  boron-  ilieon 
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Fig.  2.  THE  HIGH  CONCENTRATION  ANOMALY— COMPAR¬ 
ISON  OF  EXPERIMENTAL  DIFFUSION  PROFILE  WITH  THE 
CALCUIAiTED  PROFILE  USING  THE  DIFFUSION  COEFFI¬ 
CIENTS  DETERMINED  IN  LOW  CONCENTRATION  EXPERI¬ 
MENTS,  AFTER  NICHOLAS  [71. 

During  annealing,  the  redistribution  of  boron  is  anomalous 
At  the  beginning,  the  diffusion  rate  is  high.  However,  as  the  annealing 
proceeds,  the  diffusion  rate  decreases  to  values  obtained  in  ordinary  dif 
fusion  experiments  [2 1 .  The  relative  importance  of  this  transient  effect 
is  shown  in  Table  1,  after  Hofker  et  al  [21.  The  diffusion  coefficients 
m  the  second  and  third  columns  are  obtained  from  curve  fitting  the  exper 
imental  profiles  after  35  min  of  annealing  with  calculated  profiles  orig¬ 
inated  from  the  as-implanted  profiles.  The  values  in  the  fourth  column 
are  obtained  from  fitting  experimental  profiles  at  longer  annealing  times 
vith  calculated  profiles  originated  from  annealed  profiles.  Agreement  of 
the  values  in  the  fourth  column  with  the  diffusion  coefficients  in  column 
five,  obtained  in  ordinary  diffusions,  indicate  that,  once  the  ion  im¬ 
planted  impurities  are  on  substitutional  sites,  the  diffusion  mechanism 
is  indistinguishable  from  ordinary  diffusion. 


concentrations  under  thermal  equilibrium  conditions,  the  electrical  In¬ 
active  Boron  concentration  approaches  the  total  boron  concentration.  In 
contrast,  under  nonequilibrium  thermodynamic  conditions  such  as  at  the 
completion  of  a  boron  implant  of  a  proton  irradiation  or  before  the  com¬ 
pletion  of  an  anneal,  the  active  boron  concentration  will  differ  sub¬ 
stantially  from  the  total  boron  concentration,  as  we  will  see. 

a .  Ordinary  Diffusion 

At  low  concentrations,  the  diffusion  of  impurities  into 
near  intrinsic  silicon  substrates  obeys  Fick's  laws,  hence  the  diffused 
profile  from  a  constant  concentration  source  is  represented  by  a  comple¬ 
mentary  error  function.  At  high  concentrations,  boron  and  many  Group 
III  and  Group  V  impurities  diffuse  abnormally  fast.  This  high  concen¬ 
tration  anomaly  [3-7]  is  depicted  in  Fig.  2;  after  Nicholas  [5],  Another 
anomaly  is  the  effect  of  substrate  background  doping.  If  the  substrate 
is  doped  a-priori  with  arsenic  or  phosphorus,  the  a-posteriori  diffusion 
of  boron  is  greatly  reduced  [6,7].  These  diffusion  anomalies  can  be  rep¬ 
resented  by  the  following  equation,  after  Crowder  et  al  [8], 

D/D^  =  p/p^  (1) 

where  D  and  p  symbolize  diffusion  coefficient  and  hole  concentra t ion 
respectively,  and  the  subscript  i  denotes  intrinsic  condition.  This 
equation  can  be  derived  [8]  assuming  that  the  diffusion  of  boron  is  pro¬ 
portional  to  the  total  vacancy  concentration,  which  is  in  turn  Fermi  level 
dependent.  An  alternate  derivation  by  Anderson  and  Gibbons  [91  lias  also 
been  published  and  is,  in  fact,  consistent  with  the  annealing  model  we 
will  develop  here. 

There  are  other  anomalies,  such  as  the  enhancement  in  dif¬ 
fusion  in  mechanically  polished  substrates  and  diffusion  in  the  presence 
of  oxidation.  These  cases  require  an  alternative  mechanism  to  explain 
the  enhanced  diffusion,  for  instance,  Nicholas  [71  and  Parker  [TO]  pro¬ 
pose  that  vacancies  are  generated  from  the  nonconservative  motion  of  dis¬ 
location  . 


Next,  we  will  examine  the  experiments  on  the  lattice  location  of 
ion-implanted  boron.  These  experiments  reveal  that  at  the  initial  stages 
of  annealing  only  a  small  fraction  of  boron  is  on  substitutional  sites. 
This  fraction  gradually  increases  to  100 %  toward  the  end  of  the  anneal. 
Then,  we  review  experiments  on  the  electrical  activity  of  the  implanted 
Boron  and  show  that  this  data  is  in  good  agreement  with  the  channeling 
data  . 

It  is  evident  from  these  experiments  that  there  is  a  nonsubstitu- 
tional  fraction  of  boron  during  annealing.  Naturally,  this  inactive 
boron  need  not  be  in  a  single  type  of  defect  or  complex.  Whereas,  in 
many  cases  it  does  seem  that  a  single  type  of  defect  is  responsible  for 
the  electrical  inactivity,  in  other  cases  it  is  apparent  that  more  than 
one  form  of  inactive  boron  is  required.  To  study  this  possibility,  we 
performed  a  group  of  experiments  on  the  evolution  of  implantation  damage 
during  an  isothermal  annealing  to  elucidate  this  issue.  The  global  trans¬ 
mission  electron  microscopy  results  indicate  that,  when  the  presence  of 
more  than  one  form  of  inactive  boron  is  evident,  the  second  form  is  either 
boron  trapped  by  defects  formed  during  annealing  or  boron  precipitate,  if 
the  solid  solubility  limit  was  exceeded  by  the  boron  concentration. 

1.  Experiments  on  Diffusion 

The  study  of  diffusion  depends  heavily  on  the  determination  of 
concentration  profiles.  Whereas  some  profiling  techniques  are  based  on 
the  identification  of  the  species  chemically  present ,  others  are  based  on 
the  electrical  effect  of  the  ionized  specie  present  in  the  substrate.  In 
the  first  category,  we  find  the  use  of  radiotracers,  secondary  ion  mass 
spectrometry,  Auger,  backscatter ing ,  and  ESCA .  In  the  second  category, 
the  techniques  are  two  point  or  four  point  resistance  measurements,  CV 
profiling,  and  Hall-mobility  measurements.  In  general,  these  methods  do 
not  give  the  depth  information  unless  they  are  used  in  conjunction  with 
a  layer  removal  technique  such  as  anodic  oxidation  and  oxide  stripping 
or  sputtering  with  an  argon  ion  beam.  The  data  that  we  will  review  were 
obtained  using  SIMS  for  the  determination  of  total  boron  concentra¬ 
tion  profiles  in  the  substrate  and  Hall-mobility  measurement  to  deter¬ 
mine  the  ionized,  electrically  active  boron  profiles.  For  moderate 


temperature  under  conditions  that  include  the  most  relevant  cases  in  de- 
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vice  fabrication;  i.e.,  dose  range  of  10  to  10  ions/cm  and  annealing 
temperature  range  between  900 °C  to  1000 °C.  The  primary  input  to  the  an¬ 
nealing  model  is  the  history  of  the  substrate  summarized  by  two  distri¬ 
butions:  the  as-implanted  boron  profile  (either  computed  or  determined 

experimentally)  and  the  distribution  of  energy  deposited  into  nuclear 
processes  by  the  boron  ions.  Damage  is  a  product  of  this  energy  distri¬ 
bution.  Our  goal  is  to  use  this  basic  information  (which  is  available 
for  various  combinations  of  projectiles  and  targets  with  acceleration 
voltage  as  a  parameter)  to  predict  the  isothermal  annealing  of  boron 
implanted  into  Si  at  room  temperature  (the  most  common  application)  . 
Specifically,  we  wish  to  determine  the  evolution  of  the  electrically- 
active  boron  in  both  time  and  space.  A  realistic  annealing  model  should 
predict  ordinary  diffusion,  proton  enhanced  diffusion,  and  naturally  be 
consistent  with  the  collection  of  experiments  that  we  will  review  here¬ 
after  . 

B .  A  Summary  of  Relevant  Experimental  Data 

The  purpose  of  reviewing  relevant  experiments  is  to  establish  that, 
despite  the  complexity  of  the  actual  situation,  a  model  containing  only 
four  interacting  species  is  capable  of  explaining  the  major  features  of 
the  annealing  behavior.  These  species  are:  substitutional  boron,  boron 
vacancy  pairs,  positively  charged  vacancies,  and  immobile  boron  (B 
trapped  by  defects  or  B  precipitate).  In  what  follows,  we  will  show  that 
with  these  four  species  we  can  model  both  the  diffusive  redistribution 
of  ion  implanted  boron  and  the  electrical  activation  of  the  implanted 
boron  during  annealing. 

The  experiments  we  will  review  can  be  grouped  in  three  major  areas: 
diffusion,  electrical  activation,  and  implantation  damage.  In  the  first 
area,  we  will  examine  diffusion  anomalies,  proton  enhanced  diffusion ,  and 
the  redistribution  of  profiles  during  annealing.  A  vacancy  mechanism  is 
likely  to  account  for  the  global  features  in  these  experiments.  However, 
other  features  will  require  additional  explanations.  This  is  the  case 
for  the  abnormal  structures  in  low-temperature  and  short-time  anneals. 


shows  that  this  enhancement  in  diffusion  is  a  transient  and  that  even¬ 
tually,  as  the  electrical  activity  approaches  lOO'jt,  diffusion  will 

gradually  slow  down  and  the  ordinary  thermal  diffusion  rate  of  1.3  •< 
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10  cm  /sec  will  be  obtained. 

5.  The  Need  for  an  Annealing  Model 

The  example  in  the  preceding  section  is  from  a  common  indus¬ 
trial  application.  It  serves  the  purpose  of  showing  the  need  for  a 
realistic  annealing  model.  Such  a  model  should  provide  not  only  reas¬ 
onable  predictions  of  the  annealed  boron  profile,  but  also  a  reasonable 
prediction  of  the  evolution  of  electrical  activity  versus  time.  At 
present,  empirical  procedures  are  used  to  arrive  at  the  implantation 
conditions  and  annealing  schedule  that  will  achieve  the  desired  profiles. 
In  this  empirical  process ,  the  accurately  calculated  as-implanted  pro¬ 
files  are  only  used  as  approximations  to  the  desired  profiles.  And  the 
annealing  schedule  is-  selected  from  a  proven  set  of  conservative  sched¬ 
ules  . 

In  contrast,  given  a  suitable  annealing  model,  the  device  de¬ 
signer  can  utilize  the  calculated  as-implanted  profiles  and  from  these 
determine  the  annealed  profiles  for  one  or  more  annealing  schedules. 
Furthermore,  such  a  model  will  enable  the  entire  ion  implantation  pro¬ 
cess  to  bo  modelled  on  a  computer,  which  may  lead  to  the  eventual  devel¬ 
opment  of  more  sophisticated  applications  of  ion  implantation  in  device 
fabrication  . 

<; .  Cases  to  be  Studied 

The  principal  impurities  used  in  the  fabrication  of  silicon 
devices  are  boron  (for  p-type  regions),  phosphorus  and  arsenic  (for  n- 
type  regions).  In  this  group,  13  has  attributes  that  make  it  the  most 
widely  implanted  impurity.  For  instance,  the  implantation  of  B  for 
threshold  shifting  in  MOS  devices  and  the  implantation  of  the  boron 
doped  base  into  a  pre-prepared  substrate  with  the  E  and  C  structures 
are  possible  because  of  the  light  damage  and  the  long  range  character¬ 
istics  of  ion  implanted  boron.  We  chose  in  the  present  work  to  study 
the  annealing  behavior  of  ion-implanted  boron  into  silicon  at  room 
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thermal  equilibrium.  Therefore,  it  is  not  surprising  to  find  that  the 
redistribution  of  boron  during  annealing  cannot  be  predicted  with  a 
simple  diffusion  model  that  utilizes  diffusion  coefficients  derived 
from  ordinary  diffusion  experiments  under  thermal  equilibrium  condi¬ 
tions.  For  instance,  Fig.  1  shows  an  as-implanted  profile  for  boron 
in  silicon  together  with  this  same  profile  after  35  minutes  of  anneal¬ 
ing  at  900 °C .  Both  profiles  were  determined  experimentally  by  Hofker 
et  al  [21.  These  authors  show  that,  to  fit  the  annealed  profile  with 
the  calculated  profile  that  result  from  the  diffusion  of  the  as-im¬ 
planted  boron  distribution,  the  diffusion  coefficient  should  be  set 

2 

equal  to  1.4x10-14  cm  /sec.  This  value  is  an  order  of  magnitude  higher 
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than  1.3  x  10  cm  /sec,  the  diffusion  coefficient  determined  from  or¬ 
dinary  diffusion  experiments  extrapolated  to  900°C.  Further  annealing 
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lattice,  forming  as  a  result  a  disorder  cluster  around  the  ion  track. 
Furthermore,  since  the  process  .just  described  is  far  away  from  thermal 
equilibrium,  only  some  of  the  implanted  impurities  will  be  located  on 
substitutional  sites,  after  the  implantation.  In  general  ,  a  larger  frac¬ 
tion  will  be  found  in  interstitial  positions  and  still  other  impurity 
atoms  will  form  complexes  with  simple  lattice  defects  (lattice  vacancies 
and  interstitials).  Furthermore,  the  doping  characteristics  of  the  im¬ 
purity  may  be  overwhelmed  by  the  damage.  As  a  result,  the  as-implanted 
substrate  will  usually  exhibit  an  electrical  carrier  concentration  far 
below  that  which  would  be  predicted  simply  from  the  implanted  impurity 
concentration  profile.  In  other  words,  only  a  small  percentage  of  the 
implanted  impurities  contribute  to  electrical  carriers. 

3.  The  Anneal 

In  order  to  obtain  the  expected  electrical  properties,  it  is 
necessary  to  anneal  the  implanted  material  at  an  elevated  temperature. 
During  annealing,  some  vacancies  will  be  annihilated  by  silicon  inter¬ 
stitials,  divancies  will  anneal  out,  stress  created  by  the  damage  clust¬ 
ers  will  be  relieved  by  the  formation  of  dislocation  loops  and  vacancy 
clusters,  there  will  be  migration  of  defects,  transformation  of  com¬ 
plexes,  etc.,  and,  most  importantly,  a  relocation  of  impurities  on  sub¬ 
stitutional  sites  accompanied  by  spatial  redistribution.  After  the  an¬ 
neal,  the  substrate  will  still  retain  residual  dislocation  networks ,  but 
for  silicon  at  least,  the  mobility  and  electrical  activity  (ration  of 
electrically  active  to  the  totality  of  implanted  impurities)  are  compar¬ 
able  to  bulk-doped  samples,  and  hence  this  process  is  adequate  for  device 
fabr icat ion . 

I .  The  Spatial  Redistribution  of  Impurities  During  Annealing 

In  order  to  examine  further  the  spatial  redistribution  of  Boron 
during  the  anneal,  we  observe  that,  at  the  outset  of  the  anneal,  the  im¬ 
planted  Boron  distribution  is  in  the  presence  of  an  abnormally  high  con¬ 
centration  of  lattice  defects,  mostly  vacancies.  Consequently,  as  an¬ 
nealing  begins,  the  defects  migrate  and  transform,  and  the  redistribution 
of  boron  is  occurring  simultaneously  under  conditions  far  away  from 
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THE  ANNEALING  OF  ION  IMPLANTED  BORON  INTO  SILICON  AT  ROOM  TEMPERATURE 

A .  Introduction 

1 .  The  Problem 

Ion  implantation  is  a  technique  [1]  for  introducing  selected 
impurities  into  a  substrate  material  in  a  controlled  fashion.  It  is 
based  on  the  generation  of  a  mass  analyzed  ion  beam  followed  by  the  ac¬ 
celeration  of  these  ions  in  an  electric  field  directed  toward  the  sub¬ 
strate  material.  This  scheme  is  easily  controlled  with  electronic  in¬ 
strumentation.  As  a  result,  it  is  possible  to  achieve  repea tably  high 
purity  ion  beams  ,  to  implant  ions  to  an  accurate  dose ,  and  to  impart 
precise  kinetic  energy  to  the  ions,  thus  controlling  reproducibly  their 
stopping  points  in  the  substrate.  These  advantages  are  complemented  by 
the  availability  of  a  theoretical  model  from  which  one  can  compute  ac¬ 
curately  the  as-implanted  distribution  of  impurity  atoms.  However,  this 
latter  advantage  is  somewhat  fictitious,  because  the  implanted  ions  may 
produce  substantial  damage  in  the  (crystalline)  substrate  as  they  come 
to  rest.  To  repair  this  damage,  an  anneal  has  to  follow  the  implanta¬ 
tion.  During  this  anneal,  a  diffusive  redistribution  of  the  implanted 
Impurities  will  occur  and,  as  a  result,  the  actual  impurity  profile  may 
differ  substantially  from  the  as-implanted  profile.  The  quantitative 
prediction  of  this  redistribution  problem  is  the  principal  subject  of 
this  thesis. 

2 .  The  Implantation  of  Ions 

To  elaborate  on  the  nature  of  this  problem,  it  is  first  useful 
to  remark  that  the  distribution  of  implanted  ions  in  the  solid  substrate 
depends  on  both  the  acceleration  voltage  and  mass  relations  between  the 
implanted  ion  and  the  atoms  of  the  target.  Stopping  of  an  energetic  ion 
occurs  as  a  result  of  collisions  between  the  ion  and  the  electrons  and 
atoms  of  the  target.  If  these  collisions  are  sufficiently  energetic, 
host  atoms  will  be  displaced  from  lattice  sites,  thus  leaving  vacancies 
behind.  These  displaced  atoms  may  in  turn  displace  other  atoms  in  the 


1 


Lifetime  of  vacancies 
Positively  charged  vacancies 
Neutral  vacancy 
Space 

Spatial  increment 
Vector 

The  itJl  and  components  of  vector  y* 

th 

The  time  derivative  of  the  i  component  of 


Vll 


EA 


V+ 
f . 

l 


i.l 


K1’K2’K3’  I 
K’K'  -Kv’Kv| 
1 

P 


Pi 


SSL 

T 

T 

tdam 

Tdd 
Td  iss 


prec 
'  rel 


Electrical  activity,  fraction  of  ionized  impurities  to 
the  total  number  of  implanted  impurities 

Fermi  energy  level  in  the  band  gap 

Energy  level  of  the  positive  vacancies  in  the  band  gap 
.  th 

l  function 
Integer  index 
Integer  index 
Jacobian  matrix 

Element  in  the  Jacobian  matrix 
Constant 

Equilibrium  constant  in  the  B-BV-pair  reaction 

-5 

Boltzmann's  constant,  8.62  X  10  eV/°K 

Scaled  parameters 
Superscript  denoting  left 

_3 

Hole  concentration  in  hole/cm 

Hole  concentration  under  intrinsic  conditions 

Superscript  denoting  left 

Solid  solubility  limit 

Temperature  in  degrees  kelvin 

Lifetime  of  BV-pairs 

Time  constant  associated  with  the  release  of  vacancies 
from  the  primary  implantation  damage 

Time  constant  associated  with  the  trapping  of  boron  by 
dislocation  dipoles 

Time  constant  associated  with  the  dissolution  of  boron 
precipitates 

Time  constant  associated  with  the  precipitate  of  boron 

Time  constant  associated  with  the  release  of  boron  from 
the  annealing  dislocation  dipoles 

xi  i  i 


SYMBOLS 


B 

B~V+ 

B 

prec 

B, 


trapped 


"BI 


prec 


"E 


total 


'B 


trapped 


"BY 


"dd 


"SSL 


v  + 
eq 

D 

D1 

°B 

dbv 

D 

exp 


Ionized  boron 
Boron  vacancy  pair 
Boron  precipitate 
Immobile  boron 

Boron  trapped  by  dislocation  dipoles 

Superscript  denoting  center 

3 

Concentration  in  =/cm 

3 

Concentration  of  electrically  active  boron  in  */cm 
Concentration  of  immobile  boron 
Concentration  of  boron  precipitates 
Total  boron  concentration 

Concentration  of  boron  trapped  by  dislocation  dipoles 
Concentration  of  boron-vacancy  pairs 

Threshold  concentration  level  of  boron  for  the  formation 
of  dislocation  dipoles 

Concentration  at  the  solid  solubility  limit 

Concentration  of  positively  charged  vacancies 

Equilibrium  concentration  of  positively  charged  vacancies 

2  . 

Diffusion  coefficient  in  units  of  cm  /sec 
Diffusion  coefficient  under  intrinsic  conditions 
Diffusion  coefficient  of  electrically  active  boron 
Diffusion  coefficient  of  boron-vacancy  pairs 
Diffusion  coefficient  measured  in  an  experiment 
Diffusion  coefficient  of  vacancies 


Damage (x) 


Damage  profile  computed  from  the  distribution  of  energy 
deposited  into  atomic  processes  in  ion  implanted  silicon 


PROTON  IRRADIATION 
10  Kev,  125  nA  ,  2  hr,  550  C 


Rp 


pm 


Fig.  3.  A  COMPARISON  OF  TOTAL  BORON  PROFILKD 
DETERMINED  EXPERIMENTS  USING  SIMS  WITH  THE 
ELECTRICAL  CARRIER  PROFILE  THAT  RESULT  AFTER 
TWO  HOURS  OF  IRRADIATION  WITH  10  KeV  PROTON, 

AFTER  ANDERSON  ET  AL  ri2l. 

Hofker  et  al  [2,111.  Figure  5  displays  the  .juxtaposition  of  the  annealed 
profiles  with  the  as-implanted  boron  profiles  for  each  dose,  permitting 
the  comparison  of  the  profiles  as  the  annealing  temperature  is  varied 
(after  Hofker  et  al). 

Inspection  of  these  figures  indicates  that: 

(1)  100 J  electrical  activity  is  attained  after  35  minutes  of 

annealing  a  temperature  of  UOO°C  or  higher  and  for  im¬ 
planted  doses  of  1 0 1  ■  ’  inns /cm-  or  lower. 


(2)  In  the  10  ions/cm  cases,  the  boron  peak  concentration 
exceeds  solid  solubility.  Therefore,  there  will  be  pre¬ 
cipitation  of  boron. 

(3)  In  the  lower  dose  cases,  when  the  anneal  is  carried  out 
at  800 °C,  there  is  a  fraction  of  boron  that  is  electri¬ 
cally  inactive  and  does  not  diffuse,  hence  similar  to  a 
precipitate.  However,  the  boron  concentration  is  below 
the  solid  solubility  limit  for  these  cases. 

2 .  Experiments  on  the  Lattice  Location  of  Boron 

The  techniques  used  in  the  experiments  that  we  will  review  are 
really  indirect  measurements  that  permit  the  inference  that  a  fraction 
of  boron  is  on  substitutional  sites.  For  instance,  the  carriers  observed 
in  an  electrical  resistance  measurement  are  assumed  to  arise  from  the 
ionization  of  impurity  atoms  located  on  substitutional  sites.  However, 
substitutionality  can  also  be  measured  independently  by  the  channeling 
technique  . 

This  method  is  based  on  the  interactions  of  a  low  dose  proton 
or  Helium  ion  beam  with  the  crystal  lattice.  This  choice  of  ions  is  im¬ 
portant  because  light  mass  ions  at  low  doses  produce  negligible  damage 
to  the  crystal.  As  a  result,  the  probing  ion  will  contribute  with  lit¬ 
tle  error  to  the  measurement.  The  quantitative  measurement  is  based  on 
the  detection  of  backscattered  ions  with  a  solid  state  detector.  These 
backscattered  ions  are  the  result  of  wide  angle  collisions  of  the  light 
incident  ions  with  heavier  atoms  in  the  substrate.  During  the  measure¬ 
ment,  the  incident  beam  is  directed  along  the  channels  of  the  crystal, 
hence  the  interactions  with  the  substrate  come  primarily  from  (low  im¬ 
pact  parameter)  collisions  with  off-lattice-site  atoms  (nonsubstitutional 
atoms).  The  calculation  of  the  fraction  of  nonsubstitutional  atoms  is 
based  on  the  comparison  of  the  detector  yield  under  channelled  condition 
with  the  random  yield,  when  the  crystal  is  oriented  randomly  with  respect 
to  the  incident  beam.  The  substitutional  fraction  of  atoms  is  the  com¬ 
plementary  fraction  to  unity.  There  are  complicating  factors  that  make 
the  interpretation  more  difficult.  For  instance,  the  presence  of  lattice 
damage  and  flux  peaking  effects  [151  can  result  in  errors  that  are  sig¬ 
nificant  . 
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Channeling  Experiments 


a  . 

Figure  6a  shows  the  results  of  a  channeling  experiment  by 
North  and  Gibson  [16]  on  the  lattice  location  of  ion  implanted  boron  in 
silicon.  They  determined  that  upon  isochronal  annealing  at  a  sequence 
of  monotonically  increasing  temperatures,  the  substitutional  fraction  of 
boron  decreases  first  to  a  small  fraction  at  temperatures  between  500 °C 
and  700°C  and  then  it  increases  smoothly  to  near  90 at  higher  anneal¬ 
ing  temperatures.  Also ,  in  Fig .  6a  the  fraction  of  substitutional  boron 
can  be  compared  with  the  percentage  of  electrical  carriers  produced  by 
the  implanted  impurity  dose.  The  agreement  between  the  two  curves  is 
excellent  above  700 °C.  The  difference  at  lower  temperatures  is  attrib¬ 
uted  by  the  authors  to  the  compensation  effect  of  intrinsic  damage  pro¬ 
duced  by  the  implantation.  Although  such  an  explanation  is  reasonable 
and  transmission  electron  micrographs  on  samples  of  equivalent  prepara¬ 
tion  do  show  [17]  the  presence  of  a  high  concentration  of  vacancy  clust¬ 
ers  and  dislocation  loops,  there  is  perhaps  insufficient  attention  given 
to  the  consequences  of  dechanneling  of  the  probing  beam  by  said  defects. 
Certainly,  a  correction  of  the  data  to  include  this  effect  will  result 
in  a  lower  fraction  of  substitutional  boron  in  this  region. 

Basically,  this  experiment  indicates  that  annealing  pro¬ 
ceeds  with  the  conversion  of  nonsubst itutional ,  electrically  inactive 
boron  into  Boron  that  is  electrically  active  and  located  on  substitu¬ 
tional  sites.  A  more  detailed  inspection  of  Fig.  6a  reveals  that  between 
500°C  and  700°C  the  electrical  activity  as  well  as  the  concentration  of 
substitutional  boron  decreases  to  a  minimum.  This  phenomenon  is  known 
as  negative  or  reverse  annealing,  and  occurs  near  the  recrystal izat ion 

1  J 

temperature  of  Silicon  for  boron  implanted  in  the  dose  range  of  10±  to 
1015  ions/cm^.  Transmission  Electron  Microscopy  (TEM)  work  by  Bicknell 
[17]  indicates  that  negative  annealing  is  accompanied  by  the  formation 
of  lineated  defects  . 

The  effect  of  an  increase  in  dose  is  a  decrease  in  the 
fraction  of  substitutional  boron  and  a  lower  initial  electrical  activ¬ 
ity.  This  result  is  depicted  in  Fig.  6b  after  Fladda  et  al  (11]  and  in 
Fig.  6c  after  Seidel  and  McRae  [16]. 
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The  percentage  of  the  boron  atoms  which 
occupy  substitutional  lattice  sites  as 
a  function  of  annealing  temperature. 
Also  shown  Is  the  carrier  concentration 
expressed  as  a  percentage  of  the  maximum 
number  of  observed  carriers  determined 
from  Hall  measurements  on  a  sample  im¬ 
planted  to  a  Huence  of  2  <  10  J/em^ . 
After  North  and  Gibson  [15] . 


b.  Anneal  behavior  of  56-keV  room-tempera - 
ture  boron  Implants  In  silicon.  The 
variation  of  the  substitutional  concen¬ 
tration  with  anneal  temperature  is  shown 
for  three  different  boron  doses:  •  1.2x 
1015  ions/cm^,  A  1.6x101'*  ions/cm-,  and 
■  1.1x1014  ions/cm**.  At  a  few  selected 
temperatures,  the  substitutional  concen¬ 
tration  is  also  given  in  percent  of  the 
Implanted  dose.  After  Fladda  et  nl  [11] 


c.  Isochronal  data  In  reduced  form. 
The  ratio  of  the  free  carrier 
content  from  Hall  measurements 
to  the  dose  Is  plotted  against 
the  nnnenllng  temperature  for 
three  doses.  The  maximum  free 
carrier  content  is  nearly  #'”»al 
to  the  measured  dose  •  0. 

After  Seidel  and  MacKne  [lh]. 


Fig;.  6.  EFFECT  OF  DOSE 
ON  THE  FRACTION  OF  ION 
IMPLANTED  BORON  ON  SUB¬ 
STITUTIONAL  LOCATION 
DURING  ANNEALING. 


Transmission  Electron  Microscopy 


The  interpretation  of  transmission  electron  micrographs  is  a 
very  sophisticated  topic.  A  thumbnail  sketch  shows  two  types  of 
experimentally  observed  patterns.  A  DIFFRACTION  PATTERN  is  formed  by 
projecting  all  the  electrons  that  emerge  from  the  sample  onto  the  ob¬ 
servation  screen,  analogous  to  the  low  power  magnification  case  in  a 
light  microscope.  The  observed  pattern  consists  of  a  central  spot  and 
a  series  of  secondary  spots  and/or  rings  depending  on  whether  the  sub¬ 
strate  is  still  crystalline  or  amorphous.  These  structures  are  a  con¬ 
sequence  of  the  existence  of  or  the  lack  of  periodicity  in  the  arrange¬ 
ment  of  the  atoms  and  the  constructive  and  destructive  interference  of 
the  electron  diffraction  patterns  from  the  individual  atoms  on  the  ob¬ 
servation  screen.  The  other  type  of  pattern  is  the  DARK  or  BRIGIff  FIELD 
MICROGRAPH  which  is  obtained  by  permitting  only  those  electrons  that 
form  a  selected  dark  or  bright  area  of  the  DIFFRACTION  pattern  to  reach 
the  observation  screen.  TEM  is  an  effective  technique  for  the  study  of 
implantation  damage  and  the  annealing  behavior.  Bicknell  1 171  corre¬ 
lated  the  damage  with  electrical  activity  resulting  from  a  succession 
of  30  minute  isochronal  anneals  at  temperatures  between  dooVand  ll<>o°c. 
Naturally,  the  electrical  activity  result  follows  the  behavior  described 
previously,  namely,  a  general  increase  in  electrica 1  act  ivi  ty  interrupted 
by  reverse  annealing  between  500 °C  to  600 °c,  the  temperature  range  over 
which  amorphous  silicon  recrysta  1  i/es .  Although  no  amorphous  layer  is 
produced  in  this  experiment,  severe  migration  and  reordering  of  implan¬ 
tation  induced  defects  occur.  This  is  indicated  by  the  rapid  growth  of 
lineated  defects  in  the  micrographs.  Bicknell  interpreted  negative  an¬ 
nealing  as  precipitation  of  boron  a  tici  the  lineated  defect  as  a  product 
of  the  precipitation.  Further  annealing  at  higher  temperature  show  the 
disappearance  of  lineated  defects  and  the  emergence  of  loop  defects.  At 
even  higher  temperatures,  these  loops  grow  in  si/e  and  decrease  in  nun- 
tier.  Figure  7,  after  Bicknell,  shows  a  schematic  summary  at  the  result 
and  the  interpretation.  As  mentioned  previously,  «e  are  interested  in 
the  correlation  of  implantation  damage  with  electrical  activity  in  an 
i -otherma 1  anneal.  In  particular,  we  speculate  whether  the  lineated 
defects  will  be  present  in  isothermal  annealing  at  son  V  and  !M»o°r  for 
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Fig.  7.  SCHEMATIC  SUMMARY  OF  RE¬ 
SULTS  AND  THEIR  INTERPRETATION, 
AFTER  BICKNELL  [17]  . 


length  of  time  consistent  with  the  electrical  activation  of  boron.  Next, 
we  will  present  the  isothermal  annealing  results. 


a.  The  Isothermal  Annealing  of  Implantation  Damage 


Representative  bright  field  transmission  electron  micro¬ 
graphs  obtained  on  isothermally  annealed  samples  at  800°C  and  900°C  are 

shown  in  Figs.  8  and  9.  The  samples  were  implanted  with  70  keV  boron 
14  15  2 

to  doses  of  10  and  10  ions/cm  .  These  implantation  and  annealing 
conditions  are  identical  to  those  of  Hofker  et  al  [2,14]  and  close  to 
those  of  Seidel  and  MacRae  [18].  Consequently,  the  present  TEM  study 
[19]  can  be  correlated  with  the  electrical  activation  and  profiling  work 
by  the  previously  cited  authors.  We  can  identify  two  distinct  damage 
annealing  stages  in  the  800°C  isothermal  anneals.  The  first  stage  is 
characterized  by  the  appearance  of  a  high  concentration  of  150  to  200  A 
diameter  vacancy  dislocation  loops  inclined  the  (100)  plane,  as  shown 
in  Figs.  8b  and  9b,  at  5  minutes.  The  maximum  loop  density  occurs  at  a 
depth  of  approximately  1800  A,  in  agreement  with  calculations  of  the  lo¬ 
cation  of  the  primary  loops  is  estimated  to  be  approximately  1013  and  10^° 
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vacancies/cm  for  implantation  doses  of  10  and  10  ions/cm~,  respec¬ 
tively.  These  figures  are  about  an  order  of  magnitude  lower  than  the 
concentration  of  Frenkel  Pairs,  computed  from  the  distribution  of  energy 
deposited  by  the  implanted  boron  into  atomic  processes.  A  fraction  of 
this  difference  can  be  attributed  to  the  recombination  of  vacancies  and 
silicon  interstitials.  Since  both  calculations  in  this  comparison  are 
sub/ject  to  some  error,  we  consider  that  the  figures  are  in  reasonable 
agreement.  These  bound  vacancies  are  released  upon  annihilation  of  the 
primary  loops  between  5  and  10  minutes  of  annealing. 

The  second  stage  of  damage  annealing  appears  within  10 
minutes  of  annealing.  It  is  characterized  by  the  emergence  of  lineated 
defects  that  are  distributed  about  a  position  deeper  than  the  primary 
damage  peak  and  coincident  with  the  boron  range,  in  agreement  with  Pe- 
lous  et  al  [20] .  The  presence  of  lineated  structures  in  boron  implanted 
into  silicon  has  been  reported  by  a  number  of  authors  [21-24].  Previ¬ 
ously,  they  have  been  identified  as  rod  shaped  defects  or  boron  precip¬ 
itates  in  conflict  with  the  solid  solubility  limit.  We  interpret  the 
observed  contrast  of  these  defects  to  be  due  to  dislocation  dipoles 
[25]  and  support  Pelous's  [201  suggestion  that  the  pinning  of  boron  by 
long  range  interactions  between  boron  and  secondary  line  defects  is  a 
likely  possibility.  The  comparison  of  total  boron  concentration  and 
electrical  carrier  concentration  profiles  for  35  minutes  of  annealing 
at  800 °C  does  indicate  that  the  entire  inactive  fraction  of  boron  is 
immobile.  I f  we  associate  this  fraction  of  immobile  boron  with  the 

dislocation  dipoles*  in  Figs.  8b  and  9b  at  35  minutes,  a  simple  calcu- 

,  4 

lation  involving  defect  density  shows  that  there  are  approximately  10 

pinned  boron  atoms  per  dislocation  dipole.  Proceeding  with  this  inter¬ 
pretation,  we  reason  that  the  electrical  activation  of  boron  at  800 °c 
is  controlled  by  the  annealing  of  dislocation  dipoles.  The  900 °C  an¬ 
nealing  results  indicate  that,  as  annealing  time  increases,  the  dislo¬ 
cation  line  density  decreases  and  dissociation  of  dipoles  occur.  Large 
secondary  loops  have  been  shown  to  originate  from  dipoles,  as  illustrated 


The  residual  damage  other  than  dislocation  dipoles  in  Fig.  2b  does  not 
retain  boron,  as  inferred  by  their  presence  at  900 °C  35  min  100  '  elect. 


(circled  region  in  Fig.  9a).  Many  of  the  secondary  loops,  dislocation 
ring  structures,  and  arrays  remain  in  the  lattice  as  residual  damage 
after  35  minutes.  At  this  time,  the  electrical  activity  is  near  100', 
therefore,  this  residual  damage  does  not  trap  boron.  Presumably,  at 
900 °c,  the  two  stages  of  annealing  apparent  in  the  800 °C  cases  occurred 
at  earlier  points  in  time. 

The  experimental  evidence  on  hand  does  not  permit  a  de¬ 
tailed  explanation  on  the  exact  nature  of  how  these  dislocation  dipoles 
trap  boron  atoms.  There  are,  however,  strong  indications  that  the  for¬ 
mation  of  dislocation  dipoles  require  the  presence  of  both  a  high  con¬ 
centration  of  Boron  and  a  high  concentration  of  vacancies.  For  instance, 
the  existence  of  dislocation  dipoles  in  high  concentra t ion  diffusions 
[26 1  indicates  that  the  formation  mechanism  requires  a  boron  concentra¬ 
tion  above  a  certain  level.  We  can  also  infer  that,  in  addition  to  the 
boron,  a  large  concentration  of  vacancies  is  required  to  form  the  dis¬ 
location  dipoles.  This  deduction  follows  from  the  observation  that  the 
dislocation  dipole  density  is  much  higher  in  the  implanted  case  than  in 
the  diffused  case.  To  give  further  support  to  this  point,  we  recourse 
to  a  micrograph  in  which  various  dislocation  dipoles  are  emerging,  each 
one  from  a  vacancy  cluster.  Furthermore,  we  performed  [191  another  ir¬ 
radiation  experiment  in  which  an  annealed  ion  implanted  sample  was  sub¬ 
sequently  irradiated  with  protons  followed  by  a  short  10  minute  anneal 
at  800 .  The  transmission  electron  micrograph  of  this  sample  shows 

that,  prior  to  the  proton  bombardment,  as  described  previously,  dislo- 
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cation  dipoles  formed  and  annealed  out  at  900°C.  After  the  5  x  10 
2 

protons/cm  bombardment  at  room  temperature,  the  micrograph  shows  that 
dislocation  dipoles  were  formed  again,  presumably  when  the  vacancy  con¬ 
centration  was  greatly  increased  by  the  proton  irradiation. 

Given  the  number  of  boron  atoms  to  be  associated  with  a 

4 

dislocation  dipole  (approximately  10  boron/dd),  we  visualize  the  impur¬ 
ities  as  being  located  in  a  cloud  imbedded  in  the  strain  field  of  the 
dislocation  dipole.  Upon  annealing,  the  dislocation  dipoles  dissociate, 
releasing  the  boron  atoms.  As  mentioned  in  a  previous  section,  as  an¬ 
nealing  proceeds,  some  dislocation  dipoles  produce  secondary  defect  loops 
and  the  residual  damage  that  remains  in  the  crystal  is  ineffective  for 
trapping  boron  atoms. 


The  Annealing  Model 


C  . 

In  the  preceding  section,  we  reviewed  a  group  of  experiments  that 
encompassed  electrical  behavior  and  redistribution  of  boron  under  condi¬ 
tions  of  ordinary  diffusion,  proton  irradiation,  and  boron  introduced  in 
the  substrate  by  ion  implantation.  Before  we  develop  our  annealing  model, 
it  is  pertinent  to  highlight  some  of  the  models  proposed  to  explain  some 
of  these  experiments. 

As  mentioned  previously,  there  is  not  an  unanimous  agreement  on  the 
nature  of  the  diffusion  of  boron  in  silicon.  Details  of  this  intersti¬ 
tial-vacancy-controversy  are  well  documented  in  the  literature.  Some 
arguments  in  favor  of  an  interstitial  mechanism  can  be  found  in  an  excel¬ 
lent  review  article  by  Seeger  and  Chick  f27],  and  numerous  quantitative 
models  based  on  the  monovacancy  mechanism  extended  to  cover  different 
diffusion  anomalies  are  available  r28-30],  Basically,  a  monovacancy 
mechanism  is  a  single  stream  diffusion  model  where  the  impurity  diffu¬ 
sion  coefficient  is  linearly  dependent  upon  the  total  vacancy  concentra¬ 
tion.  Therefore,  substitution  of  the  expression  for  the  variation  of 
the  positive  vacancy  concentration  with  the  Fermi  level  in  the  bandgap 
by  Longini  and  Green  [31]  yield  an  impurity  diffusion  coefficient 
dependent  on  the  Fermi  level,  hence  on  the  impurity  concentration  and/or 
substrate  doping. 

Also  based  on  vacancies  but  in  a  different  manner,  recently  Ander¬ 
son  and  Gibbons  f9l  proposed  a  two  stream  model  for  the  diffusion  of 
boron  in  silicon.  The  two  species  are:  substitutional  boron  and  boron 
vacancy  pa ir . 

In  this  formulation,  the  diffusion  of  boron  is  a  weighted  process 
with  contributions  from  a  slow  and  a  fast  diffusing  species,  substitu¬ 
tional  boron  and  boron  vacancy  pair,  respectively.  And,  the  charged 
monovacancies  come  into  play  through  th  --  following  chemical  reaction: 

B  +  V+  =  B  V+  (2  ) 

For  instance,  an  increase  in  the  positive  vacancy  concentration  with  the 
position  of  the  Fermi  level  approaching  to  the  valence  band  edge  manifest 
a ^  a  shift  of  the  reaction  to  the  right,  hence  increasing  the  population 
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of  BY-pairs,  the  fast  diffusant.  This  increase,  in  turn,  causes  an  in¬ 
crease  in  the  overall  diffusion.  As  we  can  see,  this  approach  can  also 
predict  the  doping  anomalies,  and  it  is  possible  to  derive  from  this  for¬ 
mulation  the  basic  expression  in  Eq .  (1): 

D/Di  =  P/Pi 

This  two  stream  diffusion  model  has  also  been  shown  to  give  reasonable 
predictions  of  proton  enhanced  diffusion.  This  result  is  important  for 
two  reasons.  First,  one  would  expect  that  a  realistic  model  should  ap¬ 
ply  to  both  ordinary  diffusion  and  proton  enhanced  diffusion.  Second, 
because  1  here  are  many  similarities  between  the  PED  case  and  the  anneal¬ 
ing  problem ,  that  it  is  likely  that  the  basic  two  stream  approach  may  be 
■pplicable  to  the  annealing  behavior  of  ion-implanted  boron. 

We  -hould  also  mention  at  this  point  that  Baruch  etal  [32]  proposed 
an  ii.terst  it  ia  I  couni  erpart  of  the  two  stream  diffusion  model  for  the 
prediction  nf  proton  enhanced  diffusion.  For  ordinary  diffusion,  they 
propose  a  shift  from  the  interstitial  diffusion  mechanism  to  a  vacancy 
mocha n  i  -•••: .  Thus  far,  maybe  all  three  approaches  appear  to  be  nearly 
equivalent,  however,  when  we  examine  them  in  regard  to  the  ion-implanta¬ 
tion  related  experiments,  it  will  he  apparent  that  extending  the  basic 
two  stream  'aodel  with  boron  vacancy  pairs  to  include  positive  vacancies 
■■ay  result  in  an  annealing  model  capable  of  predicting  the  diffusive  re¬ 
distribution  of  boron  during  annealing. 

1  .  Choice  of  Diffusion  Model 

The  single  species  approach  is  incapable  of  accounting  for  both 
the  substitutional  and  nonsubs t 1 tu t iona 1  form  of  boron.  In  the  two 
stream  diffusion  approach,  the  BV-pair  or  boron  interstitial  are  natu¬ 
ral  lv  a  •  -oc iated  with  the  nonsubs t i tut iona 1  boron.  However,  internal 
frieti  n  I  1  and  phot  •  ic.indiirt  iv  i  ty  experiments  |\'H  1  have  indicated  that 
boron  interstitials  anneal  out  near  Boo  C:  therefore,  it  is  unlikely  that 
they  ire  nn--ent  at  eon  c.  Rerent  lv,  Watkins  [3  a  1  reported  tile  experi¬ 
mental  i”,  idi’se  of  .m  Electron  Paramagnetic  Resonance  (KPR )  peak  that  ho 
tentative]  v  identified  related  t<>  BY -pa  i  v~  .  However,  in  his  expert  - 
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THE  ANNEALING  CASES 
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(a)  The  species  are:  substitutional  boron,  electrically 
active  and  a  slow  diffusant;  BV-pn ir ,  electrically 
inactive  and  a  fast  diffusant,  and  positive  vacancies. 

(b )  A  reaction  exists  between  B,  BV-pairs.  and  positive 
vacancies  .  The  kinetics  of  this  reaction  governs  the 
time  evolution  of  the  populations  of  active  boron  and 
BV-pairs.  Therefore,  it  also  governs  the  electrical 
activation  and  indirectly  the  reci  istribution  . 

fi-l  The  diffusion  of  boron  is  a  weighted  process  with  con¬ 
tributions  of  BV-pairs,  the  fast  diffusant,  and  sub- 
situtional  boron,  the  slow  diffusant. 

(d  )  I'nder  thermal  ei|U  i  1  ibr  ium  ,  the  concentration  of  posi¬ 
tive  vacancies  is  a  function  of  the  temperature  and 
the  Fermi  level  in  the  bandgap. 


H  .  11  igh  Dose  Ca  sc 

The  high  do-e  cases  are  represented  in  the  left-hand  column  in 

1  c> 

Table  A.  The  canonical  case  representing  this  group  is  the  I11  ions-cm 

dose  implant  ,  anneal od  j  t  boo  V .  In  this  case,  the  peak  boron  coneentra 
do  ,'i 

tion  is  m  a  10  atoms  G-r:  and  the  solid  solubility  limit  at  DOO  (’  is  i. 
20  3 

lo  a  toms  cm  ;  therefore,  we  expect  the  boron  atoms  fo  precipitate.  The 


(c)  Under  therr.nl  equilibrium,  because  the  lifetime  of  va¬ 
cancies  is  short  compared  to  other  time  constants  in  the 
diffusion  process,  it  is  possible  to  approximate  that  the 
vacancy  concentration  is  the  thermal  equilibrium  concen¬ 
tration  . 

(d )  During  the  annealing  of  room  temperature  boron  implants, 
it  is  possible  to  model  the  vacancy  perturbation  arising 
from  the  implantation  as  a  vacancy  pulse  superimposed  on 
a  background  vacancy  concentration  at  equilibrium  with 
the  substrate. 

(e)  In  many  instances,  the  effect  of  the  vacancy  pulse  (aris¬ 
ing  from  the  annealing  of  vacancy  clusters  and  vacancy 
loops)  on  the  overall  annealing  is  not  significant.  This 
effect  is  only  apparent  in  short  anneals  at  high  anneal¬ 
ing  temperature. 

(f)  This  description  of  vacancies  does  not  apply  to  the  an¬ 
nealing  of  boron  that  is  implanted  at  liquid  nitrogen 
temperature  or  the  case  of  the  implantation  of  a  heavy 
ion  . 


"> .  Summary  of  the  Typical  Annealing  Cases 

Within  the  group  represented  by  this  canonical  case,  annealing 
a  ill  proceed  as  described  with  the  following  trends  as  temperature  and 
implantation  dose  are  varied.  An  increase  in  annealing  temperature  will 
result  in  a  decrease  in  the  reaction  time  constant.  Consequently,  an¬ 
nealing  will  be  faster  (a  shorter  transient  is  equivalent  to  a  faster 
attainment  of  equilibrium  conditions).  An  increase  in  implantation  dose 
produces  proportionally  higher  boron  concentrations  but  lower  initial 
electrical  activity  as  observed  by  Fladdn  et  al  [111  and  Seidel  and  Mac 
hue  rtHl .  Therefore,  annealing  time  will  be  longer  than  for  a  lower 
l"'!'  ease  annealed  at  the  same  temperature.  This  situation  holds  true 
until  the  increased  dose  causes  a  peak  boron  concentration  to  exceed  the 
-"lid  solubility  limit  at  the  annealing  temperature.  Under  these  cir¬ 
cumstances,  annealing  will  no  longer  obey  the  description  n i von  above, 
m~tead  it  .ill  follow  the  global  features  of  high  dose  cases  ropre- 
t-nted  in  the  left-hand  column  in  ^able  2. 

We  conclude  the  discussion  of  the  first  canonical  case  which 
i -  representative  of  the  cases  in  the  lower  left-hand  corner  < 
hv  sunnari/ing  that: 
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if  Table  A 


In  this  expression,  Damage(x)  is  the  vacancy  distribution  computed 

from  the  distribution  of  energy  deposited  into  nuclear  processes,  k  is 

a  constant  adjusted  to  produce  the  total  number  of  vacancies  calculated 

in  Fig.  11,  and  x,  in  the  exponential  factor  determines  the  duration 
dam 

of  the  generation  of  vacancies. 

In  this  section,  we  modelled  three  processes  that  take  place 
during  the  annealing  of  room  temperature  implantation  of  boron  into 
silicon.  These  processes  are:  the  release  of  vacancies  from  vacancy 
clusters  and  vacancy  loops,  the  absorption  of  vacancies  by  the  bulk  and 
the  surface  to  restore  the  equilibrium  vacancy  concentration  and  the 
spontaneous  breakup  of  BV-pairs.  During  the  anneal,  the  great  majority 
of  vacancies  generated  are  absorbed  by  the  substrate  and  the  smaller 
remaining  fraction  may  increase  the  population  of  BV-pairs  only  if  the 
time  constant  associated  with  the  generation  of  vacancies  is  comparable 
to  the  lifetime  of  BV-pairs.  At  annealing  temperatures  below  1000°C, 
the  lifetime  of  BV-pairs  is  long  compared  to  the  duration  of  the  gener¬ 
ation  of  vacancies.  Consequently,  the  boron-BV-pa ir  reaction  finds 
essentially  the  concentration  of  positive  vacancies  at  equilibrium.  At 
1000 °c  and  above,  the  lifetime  of  BV-pairs  is  short;  therefore,  there 
is  some  overlap  between  the  generation  of  vacancies  and  the  breaking  up 
of  BV-pairs.  For  this  reason,  the  shape  of  the  annealed  profile  is  mod¬ 
ified  by  the  localized  enhancement  in  diffusion  caused  by  the  vacancies 
that  arise  from  the  annealing  of  vacancy  clusters  and  loops.  This  effect 
is  only  observable  at  short  annealing  times.  At  longer  annealing  times, 
it  is  masked  by  ordinary  diffusion  which  is  more  important  at  inon°c  and 
above.  It  is  appropriate  to  emphasize  that  this  description  pertains  to 
the  annealing  of  a  silicon  substrate  implanted  with  boron  ions  at  room 
temperature.  The  behavior  of  vacancies  under  implantation  conditions 
that  produce  an  amorphous  layer  will  certainly  be  different. 

The  major  results  in  this  section  can  be  summarized  as: 

(a)  Fnder  thermal  equilibrium,  the  concentration  of  neutral 
vacancies  is  a  function  of  temperature  only  [Bfil. 

(b)  Under  thermal  equilibrium,  the  concentration  of  positive 
vacancies  is  a  function  of  temperature  and  the  Fermi- 
level  in  the  bandgap  T3B1. 
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Fig.  11.  THE  CALCULATED  VACANCY  CON¬ 
CENTRATION  PROFILE  PRODUCED  BY  THE 
IMPLANTATION  OF  BORON. 


above  the  equilibrium  vacancy  concentration.  Consequently,  the  substrate 
is  severely  stressed  and  it  is  incapable  of  sustaining  such  an  extraor¬ 
dinary  excess  vacancy  concentration.  As  a  result,  the  vacancies  coalesce 
in  vacancy  clusters  and  vacancy  loops  to  relieve  the  stress  in  the  lat¬ 
tice.  Transmission  electron  micrographs  reveal  that  bound  vacancies  will 
be  released  from  the  vacancy  clusters  and  vacancy  loops  as  these  defects 
anneal  out  and  that  the  background  vacancy  concentration  will  reach 
thermal  equilibrium  very  rapidly.  This  result  is  in  agreement  with  the 
backscattering  work  by  Westmoreland  et  al  [391.  The  process  we  iust  de¬ 
scribed  can  be  modelled  with  the  vacancy  generation  pulse  in  Eq .  (Id), 
that  we  propose. 

Vacancy  Generation  Rate  =  k  Damage(x)  exp(-t/T  )  (10) 

DAM 
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E  =  energy  level  of  positive  vacancies,  0.35  eV  above  the  valence 
^  band  edge 

Ep,  =  Fermi  level 

K  =  Boltzmann  constant 

T  =  temperature  in  degrees  Kelvin 

The  Fermi  level  is,  in  turn,  dependent  on  the  concentration  of  active 
boron  through  the  Fermi-Dirac  statistics  .  It  is  apparent  now  that  the 
Fermi  level  dependence  of  the  concentration  of  positive  vacancies  is  the 
mechanism  responsible  for  the  high-concentration  and  substrate-doping 
anomalies  observed  in  ordinary  diffusions  under  thermal  equilibrium  con¬ 
ditions  . 

Under  nonequilibrium  conditions,  if  the  perturbation  is  small, 
the  bulk  and  the  surface  of  the  semiconductor  will  absorb  or  generate 
vacancies  to  restore  the  equilibrium  vacancy  concentration.  We  model  the 
rate  for  the  attainment  of  equilibrium  using  the  departure  from  equilib¬ 
rium  as  the  driving  force  and  a  time  constant  t  ,  see  Eq .  (9): 


Rate  for  the  Attainment  of  Equilibrium  = 


where  is  the  vacancy  lifetime.  The  values  that  Ty  take  are  small 

compared  to  other  time  constants  associated  with  diffusion.  Consequently, 
in  those  cases  where  the  vacancies  are  only  slightly  perturbed,  it  is 
possible  to  approximate  that  the  vacancy  concentration  is  the  equilibrium 
concentration  . 

Naturally,  when  the  perturbation  is  large,  the  aforementioned 
approximation  may  no  longer  be  adequate.  For  instance,  Fig.  11  shows  the 
vacancy  concentration  produced  by  the  implantation  of  boron  ions  at  70 
keV  into  silicon.  This  profile  is  computed  from  the  number  of  Frenkel- 
pairs  produced  by  the  implantation;  which  is,  in  turn,  calculated  from 
the  distribution  of  energy  deposited  into  atomic  processes  by  the  ener¬ 
getic  ions  r 3 S 1  and  by  assuming  that  12  eV  are  required  to  produce  a 
Frenkel -pa ir .  The  calculated  vacancy  profile  is  many  orders  of  magnitude 


Reaction  Rate 


(5) 


CBV  -  koCBCV 


This  term  will  be  used  in  the  mathematical  formulation  of  the  annealing 
model.  The  time  constant  r  determines  the  rate  of  the  BV-pair-substi- 
tutional  boron  reaction.  Therefore,  indirectly  it  also  determines  the 
rate  of  the  anneal  and  the  diffusive  redistribution. 


4 .  The  Positively  Charged  Vacancies 

The  presence  of  positive  vacancies  in  this  model  has  important 
implications.  For  instance,  rearranging  Eq.  (4),  we  obtain: 


BV  ,  „ 

CB  -  koV 


(6) 


In  this  expression,  we  can  appreciate  that  the  concentration  of  positive 
vacancies  determines  at  equilibrium  the  relative  concentrations  of  active 
boron  and  BV-pairs,  hence  the  overall  diffusion. 

Vacancies  exist  in  the  silicon  lattice  in  various  charged 
states.  The  neutral  vacancy  is  a  nonionized  flaw;  therefore,  according 
to  the  analysis  by  Shockley  and  Moll  [36],  its  equilibrium  concentration 
is  a  function  of  the  temperature  only.  We  use  the  expression  in  Eq .  (7), 
after  Seidel  and  MacRae  [37] 

CyQ  =  5  1025  exp (-1.7  eV/KT)  (7) 

eq 

The  positive  vacancy  is,  on  the  other  hand,  an  ionized  flaw.  Conse¬ 
quently,  its  equilibrium  concentration  is  a  function  of  both  temperature 
and  the  Fermi  level  in  the  bandgap: 


C  (T) 
V° 


where 


C  _  =  concentration  of  neutral  vacancies 
V° 


[Eq  . 


(7)1 


(») 
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concentration  above  the  equilibrium  level.  Fladda  et  al  Til]  observed 
that  a  Riven  dose  of  protons  removes  a  fixed  amount  of  boron  from  sub¬ 
stitutional  sites,  rendering  the  boron  electrically  inactive.  Further¬ 
more,  they  observed  that  doubling  the  proton  dose  will,  in  turn,  double 
the  number  of  boron  atoms  removed  from  substitutional  sites.  This  cor¬ 
respondence  between  the  number  of  boron  atoms  rendered  inactive  with  the 
proton  dose  or  indirectly  with  the  number  of  vacancies  produced,  leads 
us  to  infer  the  existence  of  a  chemical  react  ion  between  a  substitutional 
boron  and  a  vacancy.  Since  BV-pairs  are  the  inactive  boron,  it  is  ap¬ 
parent  that  a  substitutional  boron  reacts  with  a  (positive)  lattice  va¬ 
cancy  to  form  a  BV-pair.  Because  the  BV-pair  is  electrically  neutral 
and  the  substitutional  boron  is  ionized  negatively,  the  presence  of  the 
positive  vacancy  is  necessary  in  order  to  preserve  the  conservation  of 
charge  in  the  proposed  reaction: 


B_  +  V+  ^  B~V+  (3) 

Let  us  elaborate  on  the  implications  of  this  reaction.  We  just  saw  that 
a  proton  implant  drives  the  vacancy  concentration  above  the  equilibrium 
level  forcing  the  reaction  in  Eq .  (3)  to  the  right.  And,  we  recall  that 
during  annealing  the  high  BV-pair  concentration  forces  the  reaction  in 
Eq .  (3)  to  the  left.  When  the  reactions  in  opposing  directions  cancel 
each  other,  the  reaction  is  at  equilibrium.  Under  this  special  circum¬ 
stance,  a  simple  algebraic  expression  exists  relating  the  concentrations 
of  reactants  and  product: 


BV 


=  k  C  C  , 
o  B  V+ 


(4) 


where  k  is  the  equation  constant,  a  function  of  temperature.  Whereas 
o 

the  equality  of  the  members  in  this  equation  is  indicative  of  equilib¬ 
rium,  the  inequality  is  a  measure  of  the  departure  from  thermodynamic 
equilibrium.  Consequently,  the  kinetics  of  this  reaction  can  be  mod¬ 
elled  by  the  driving  force  represented  by  the  departure  from  equilibrium 
and  a  time  constant  t,  as  indicated  n  Eq .  (5): 
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of  BV-pairs,  and  enhanced  diffusion,  see  Fig.  10b.  On  the  other  hand, 
at  the  completion  of  the  anneal,  we  find  a  thermal  equilibrium  situation 
characterized  by  high  electrical  activity,  a  large  population  of  substi¬ 
tutional  boron,  and  ordinary  diffusion.  We  insert  these  dichotomous  at¬ 
tributes  in  the  figure  as  mnemonics  for  our  next  argument.  By  inspection 
of  Fig.  10b,  we  can  interpret  annealing  as  a  transient  process  going  from 
a  nonequilibrium  state  toward  an  equilibrium  state  in  which  a  large  ini¬ 
tial  population  of  BV-pairs  transform  into  electrically  active  substitu¬ 
tional  boron.  The  comparison  of  the  diffusive  attributes  with  the  popu¬ 
lation  of  BV-pairs  and  substitutional  boron  at  the  outset  and  end  of  the 
anneal  lead  us  to  infer  that  BV'-pairs  are  fast  diffusants  and  substitu¬ 
tional  boron  are  slow  diffusants.  With  this  interpretation,  the  diffu¬ 
sive  behavior  during  annealing  becomes  a  weighted  process  with  contribu¬ 
tions  from  a  fast  and  electrically  inactive  species,  the  BV-pair,  and  a 
slow  and  electrically  active  fraction,  the  substitutional  boron.  Hence, 
a  large  population  of  BV-pairs  lead  to  enhanced  diffusion,  whereas  a 
large  population  of  substitutional  boron  produces  ordinary  diffusion. 

In  summary,  annealing  is  governed  by  the  driving  forces  that 
arise  from  a  nonequilibrium  state.  These  driving  forces  produce  the 
spontaneous  conversion  of  BV-pairs  into  substitutional  boron.  And,  the 
electrical  activation  and  diffusive  behavior  are  consequences  of  the  di¬ 
chotomous  attributes  of  BV-pairs  and  substitutional  boron  with  respect 
to  electrical  activity  and  diffusion. 

3 .  The  'Boron-Boron  Vacancy  Pair'  Reaction 

In  the  preceding  section,  we  concluded  that  the  spontaneous 
conversion  of  BV-pairs  into  substitutional  boron  governed  the  annealing 
behavior.  Now,  we  will  examine  this  conversion  problem  in  greater  de¬ 
tail.  i  le  transformation  of  BV-pairs  into  substitutional  boron  comes 
gradually  to  an  end  as  the  equilibrium  conditions  are  attained.  We  can 
now  address  the  question  whether  it  is  possible  to  perturb  the  equilib¬ 
rium  condition  and  force  the  transformation  in  the  reverse  direction? 
The  answer  is  yes.  It  is  possible  to  perturb  the  equ i 1 ibr ium  by  implant¬ 
ing  protons.  Protons  are  light  ions.  Therefore,  the  damage  they  produce 
is  also  light,  consisting  mainly  of  simple  defects  that  raise  the  vacancy 


this  particular  case,  it  does  appear  that  all  the  inactive  boron  is  in 
the  form  of  boron-vacancy  pairs.  During  annealing,  as  we  can  appreciate 
from  observing  the  monotonic  nature  of  the  electrical  activity  versus 
time  curve,  it  is  possible  to  relocate  the  nonactive  boron  on  substitu¬ 
tional  sites,  hence  increasing  the  electrical  activity  to  near  \00%  in 
35  min.  During  annealing,  a  succession  of  events  take  place:  there  is 
annihilation  of  vacancies  by  interstitials,  annealing  of  divacancies, 
formation  of  vacancy  clusters  and  vacancy  dislocation  loops  as  a  means 
of  releasing  stress  created  by  the  vacancy  clusters  in  the  crystal  lat¬ 
tice,  and,  most  importantly,  there  is  relocation  of  inactive  boron  on 
substitutional  sites  accompanied  by  redistribution.  This  conversion  of 
electrically  inactive  boron  to  electrically  active  boron  is  represented 
in  Fig.  10b  by  the  arrow  crossing  the  electrical  activity  curve  which  is 
the  boundary  between  the  two  forms  of  boron.  As  the  annealing  proceeds, 
the  fraction  of  electrically  active  boron  increases  gradually  while  the 
impurity  atoms  diffuse  in  the  crystalline  substrate.  Finally,  at  the  end 
of  the  35  min,  we  find  that  the  implanted  distribution  has  displaced  sig¬ 
nificantly  and  it  is  now  represented  by  the  deeper  boron  profile  in  Fig. 
10a.  An  attempt  to  predict  this  diffusive  redistribution  would  be  to 

take  the  diffusion  coefficient  of  boron  in  silicon  at  the  annealing  tem- 
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perature  of  900°C,  which  is  1.4  10  cm  /sec,  and  use  this  number  to 
calculate  where  the  implanted  profile  will  diffuse  to  after  35  minutes 
have  elapsed .  This  calculation  shows  that  the  profile  would  move  only 
very  slightly.  Hofker  et  al  [2]  show  that,  in  order  to  fit  the  experi¬ 
mental  profile  obtained  after  the  annealing,  the  diffusion  coefficient 

-14  2 

in  the  aforementioned  calculation  has  to  be  increased  to  1.7  10  cm  / 
sec,  which  is  one  order  of  magnitude  higher  than  the  ordinary  diffusion 
rate.  These  calculations  clearly  indicate  the  presence  of  enhanced  dif- 
fusion  during  annealing.  Further  annealing  indicates  that  the  enhanced 
diffusion  is  only  transitory  because,  as  the  electrical  activity  ap¬ 
proaches  to  lOO'i,  the  enhanced  diffusion  reduces  to  ordinary  diffusion. 

2  .  The  Spontaneous  Conversion  of  BV-Pairs 

At  the  outset  of  the  anneal  ,  we  encounter  a  nonequi librium 
situation  characterized  by  low  electrical  activity,  a  large  population 
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Fig.  10.  THE  ELECTRICAL  ACTIVATION  AND  REDISTRIBUTION 
BEHAVIOR  OF  BORON  IN  TYPICAL  ANNEALING  CASES . 


i 

The  process  we  .just  described  is  far  away  from  thermal  equi¬ 
librium.  Therefore,  only  a  small  fraction  of  boron  atoms  will  be  located 
on  substitutional  sites,  a  larger  fraction  will  be  on  interstitial  sites, 
^  and  other  boron  atoms  will  form  complexes  with  simple  defects,  i.e.,  va¬ 

cancies  and  silicon  interstitials.  Consequently,  the  electrical  activity 
at  the  completion  of  the  implant  is  far  below  what  one  would  expect  from 
the  implanted  dose.  This  situation  is  depicted  in  Fig.  10b  at  the  origin 
of  the  time  axis,  in  the  horizontal  direction.  The  electrical  activity 
is  represented  on  the  vertical  axis.  At  any  instant  in  time,  the  frac¬ 
tion  of  boron  that  contributes  to  electrical  carriers  is  represented  by 
the  ordinate  of  the  electrical  activity  curve.  On  the  other  hand,  the 
nonelectrically  active  boron  is  represented  above  the  curve  by  the  dis¬ 
tance  between  the  100  '  activity  line  and  the  ordinate  of  the  curve. 
Therefore,  we  can  also  interpret  Fig.  10b  as  a  bar  graph  describing  the 
relative  concentrations  of  active  and  inactive  boron.  In  general,  the 
nonelectrically-active  boron  need  not  be  in  a  single  form;  however,  in 
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produced  by  the  10  ions/cm  doses  exceed  the  solid  solubility  limit  at 
all  the  annealing  temperatures  in  the  table.  Therefore,  boron  will  pre¬ 
cipitate.  And,  the  precipitated  phase  is  characteristic  of  the  high  dose 
implantation  cases.  The  remaining  cases  are  represented  in  the  upper 
left-hand  corner  of  the  table  and  they  are  characterized  by  the  trapping 
of  boron  atoms  by  dislocation  dipoles  that  form  during  the  early  stages 
of  annealing.  At  800°C,  these  dislocation  dipoles  do  not  anneal  out 
completely  in  35  minutes.  Therefore,  they  are  incapable  of  releasing 
the  boron  atoms  trapped  in  the  associated  strain  fields. 

These  three  groups  have  distinct  characteristics  that  wil]  manifest 
in  the  electrical  activation  as  well  as  in  the  redistribution  of  the  im¬ 
planted  boron.  In  the  discussions  that  will  follow,  we  will  select  can¬ 
onical  cases  representative  of  each  of  the  aforementioned  groups.  Our 
exposition  will  concentrate  on  each  of  the  three  canonical  cases  and 
on  the  global  features  in  the  anneal  when  implantation  dose  and  anneal¬ 
ing  temperature  are  varied. 


1 .  The  Typical  Annealing  Cases 

14  2 

We  first  examine  the  10  ions/cm  dose,  70  keV  boron  implant 
annealed  at  900 °C  for  35  minutes.  This  is  a  relatively  simple  case.  In 
the  insert  in  Fig.  10a,  we  depict  the  silicon  surface  coincident  with 
the  origin  of  the  horizontal  axis  representing  depth  in  tenth  of  microns 
into  the  substrate  material.  The  70  keV  boron  beam  is  directed  toward 
the  surface  from  left  to  right.  As  these  energetic  ions  penetrate  the 
target,  their  stopping  in  the  substrate  is  produced  by  collisions  with 
electrons  and  host  atoms.  In  this  process,  some  collisions  are  suffi¬ 
ciently  energetic  that  host  atoms  are  displaced  from  their  lattice  sites 
leaving  vacancies  behind;  the  displaced  atoms  may  in  turn  displace  other 
atoms.  This  process  leaves  a  disorder  cluster  around  the  ion  track.  The 
insert  in  Fig.  10a  represents  one  disorder  cluster.  Since  this  collision 
process  is  random  in  nature,  not  all  the  implanted  boron  ions  will  stop 
at  the  same  depth.  Instead,  the  ions  will  come  to  rest  at  different 
depths  with  different  probabilities.  Thus,  after  the  implant,  the  dis¬ 
tribution  of  implanted  ions  is  represented  by  the  impurity  concentration 

18  2 

profile  shown  in  Fig.  10a.  The  peak  concentration  is  8  v  io  atoms/cm 
14  2 

for  the  10  ions/cm“  dose  implant  in  this  example. 
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Table  2 


ISOCHRONAL  ANNEALING  RESULTS 
Annealing  Time  =  35  min 


Dose:  ions/cm 


When  the  annealing  is  performed  isochronally  for  35  min,  in  this  range 
of  temperatures  the  electrical  activity  (at  the  completion  of  the  an¬ 
neal)  will  vary  from  a  small  fraction  to  near  100^,  depending  on  the 
particular  implantation  dose  and  annealing  temperature.  In  other  words, 
the  ratio  of  implanted  boron  that  contribute  electrical  carriers  (holes) 
to  the  totality  of  implanted  boron  ions  is  a  function  of  implantation 
dose,  annealing  temperature,  and  time.  The  understanding  of  these  re¬ 
lationships  is  one  of  our  principal  objectives. 

This  collection  of  implantation  doses  and  annealing  temperatures 
illustrated  in  Table  2  covers  a  very  wide  range  of  conditions.  For  in¬ 
stance,  most  applications  of  ion  implantation  in  device  fabrication  are 
performed  under  conditions  represented  in  the  lower  left-hand  corner  of 
this  table.  These  cases  are  very  relevant  because  the  electrical  activ 
ities  approach  to  near  100  '  after  35  minutes  of  ar  ing .  Among  the 
remaining  cases  in  this  table,  we  can  distinguish  two  groups.  The  first 
group  is  composed  of  high  dose  implants  represented  by  the  right-hand 
column  of  the  table.  In  these  cases,  the  peak  impurity  concentration 


in  order  to  predict  the  annealing  behavior  under  all  the  annealing  tem¬ 
perature  and  implantation  dose  ranges  that  we  are  interested  in. 


2 .  The  Immobile  and  Nonsubstitutional  Boron 

A  comparison  of  the  electrical  carrier  concentration  profiles 
and  the  total  boron  profiles  for  various  implantation  doses  and  anneal¬ 
ing  temperatures  is  depicted  in  Fig.  4.  It  is  apparent  from  the  inspec¬ 
tion  of  these  cases  that,  in  the  high  dose  and/or  low  annealing  tempera¬ 
ture  cases,  there  is  a  fraction  of  electrically  active  boron  that  is 
IMMOBILE.  This  new  fraction  is  present  near  the  peak  of  the  boron 
profile.  In  the  high  dose,  10^  ions/cm^  cases,  the  boron  peak  exceeds 
the  solid  solubility  limit;  therefore,  the  immobile  boron  is  a  precipi¬ 
tate.  However,  in  the  moderate  dose  cases,  the  peak  boron  concentration 
does  not  exceed  the  solid  solubility  limit;  hence,  this  immobile  boron 
is  not  a  precipitate  in  the  conventional  sense.  We  propose  that  this 
immobile  boron  consists  of  boron  atoms  trapped  by  dislocation  dipoles 
that  form  during  the  annealing. 


D .  Development  of  the  Model 

In  the  preceding  sections,  we  reviewed  a  group  of  experiments  re¬ 
lated  to  ion-implantation  and  our  interpretation  of  these  experiments. 

In  this  section,  we  will  use  the  four  species  inferred  previously  (sub¬ 
stitutional  boron,  boron-vacancy  pairs,  positive  vacancies  ,  and  immobile 
boron)  to  formulate  a  model  capable  of  describing  the  overall  annealing 
behavior  of  boron  implanted  into  silicon  at  room  temperature.  Specifi¬ 
cally,  we  wish  to  explain  the  redistribution  and  the  electrical  activa¬ 
tion  of  the  implanted  boron.  In  the  course  of  this  discussion,  the  key 
attributes  in  the  annealing  will  be  apparent  and  the  global  features 
that  result,  when  the  implantation  dose  and  anneal  temperature  are  var¬ 
ied,  will  also  follow  naturally. 

The  cases  to  be  studied  are  shown  schematically  in  Table  2.  In  the 
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horizontal  direction,  we  show  boron  implantation  doses  varying  from  10 
lfi  2 

to  10  ions/cm  .  For  a  70  keV  boron  implant,  this  dose  range  will  lead 

19  21  ,3 

to  peak  concentrations  between  10  to  10  atoms/cm  .  In  the  vertical 
direction,  we  show  annealing  temperature  ranging  between  800 °c  to  1000 °c . 


We  u  e  faced  with  a  difficult  decision.  Until  more  evidence 
on  the  stability  and  positive  identification  of  the  BV-pair  is  availa¬ 
ble,  we  proceed  with  the  examination  of  the  boron-BV-pa ir  two  stream 
diffusion  models  extended  to  include  the  diffusion  of  vacancies  .  This 
extension  is  necessary  because,  in  general,  the  migration  of  vacancies 
during  annealing  cannot  be  safely  neglected.  In  this  three-stream 
diffusion  model,  the  electrical  activity  becomes  the  ratio  of  substitu¬ 
tional  boron  to  the  sum  of  BV-pairs  and  B  substitutional  .  The  stoichi¬ 
ometry  of  the  reaction,  namely,  one  substitutional  to  boron  and  one 
vacancy  to  form  one  BV-pair  can  correspond  to  the  fixed  proportion  of 
boron  removed  from  substitutional  sites  per  unit  dose  of  proton  irradi¬ 
ation  in  the  experiments  by  Fladda  et  al  [11],  Furthermore,  let  us  ex¬ 
amine  qualitatively  if  the  diffusion  attributes  of  the  substitutional 
boron  and  BV-pair  are  consistent  with  the  red istribution  of  ion-implanted 
boron  during  annealing.  The  small  fraction  of  substitutional  boron  at 
the  outset  of  the  anneal  implies  a  large  fraction  of  BV-pairs,  the  fast 
diffusant.  Annealing  proceeds  then  with  the  conversion  of  nonsubstitu- 
tional  boron  or  BV-pairs  into  substitutional  boron  as  indicated  by  the 
experiments  on  the  lattice  location  of  boron.  Therefore,  the  redistri¬ 
bution  during  annealing  is  characterized  by  an  initial  fast  diffusion 
caused  by  a  large  fraction  of  BV-pairs  reducing  gradually  to  a  slow  dif¬ 
fusion  when  the  population  of  substitutional  boron  predominates.  This 
prediction  is  in  agreement  with  the  anomalous  transient  diffusion  ob¬ 
served  in  annealing  experiments.  Qualitatively,  there  is  an  obvious 
relation  between  the  overall  diffusion  coefficient  and  the  electrical 
activity.  Furthermore ,  if  the  chemical  reaction  in  Eq .  (2)  is  the  dom¬ 
inant  annealing  mechanism,  then  it  is  possible  to  relate  the  lifetime 
of  the  BV-pairs  to  the  annealing  time  at  a  particular  temperature.  This 
is  the  kind  of  reasoning  that  we  use  to  determine  the  values  of  the  pa¬ 
rameters  in  the  model;  some  of  them  not  determined  previously,  others  we 
will  obtain  by  alternative  means. 

In  summary,  the  three  stream  diffusion  model  appears  to  contain 
the  global  features  for  explaining  the  annealing  behavior  of  ion  im¬ 
planted  boron  into  silicon.  We  shall  see  later  that  it  will  be  necessary 
to  extend  this  model  to  include  other  forms  of  nonsubstitutional  boron 
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experimental  annealing  result  shows  in  Fig.  111a  that  there  is  diffusion 
of  the  tail  in  the  impurity  concentration  profile,  but  the  central  por¬ 
tion  near  tile  peak  remained  unchanged.  Furthermore,  in  the  same  figure 
we  can  appreciate  from  the  inspection  of  the  annealed  electrical  carrier 
profile  that  the  immobile  boron  is  also  electrically  inactive.  At  the 
outlet  of  the  anneal,  especially  at  this  high  dose,  the  electrical  ac¬ 
tivity  is  quite  low  0-12?);  therefore,  nearly  all  the  implanted  boron 
is  electrically  inactive.  From  these  facts,  we  conclude  that  the  com¬ 
position  of  tiie  inactive  boron  is  mostly  boron  precipitate  and  a  much 
smaller  fraction  of  BV-pairs.  The  experimental  evidence  on  hand  does 
not  indicate  whether  the  boron  precipitate  is  formed  directly  from  the 
implanted  boron  once  the  solid  solubility  limit  is  exceeded  or  whether 
BV-pairs  are  first  formed  and  then  followed  by  precipitation.  We  choose 
tin-  second  possibility  for  the  sake  of  convenience  in  our  mathematical 
model.  This  arbitrary  choice  does  not  modify  the  end  result  because 
precipitation  takes  place  very  rapidly  regardless  of  the  boron  source, 
t'onse  ,’iently,  in  this  interpretation,  annealing  of  this  canonical  case 
••■ill  proceed  as  follows.  At  the  outset,  the  low  activity  is  due  to  a 
large  fraction  of  BV-pairs,  from  which  the  fraction  above  the  solid 
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solubility  limit  precipitates  very  rapidly.  After  a  short  time,  we  find 
the  situation  depicted  in  the  insert  in  Fig.  12b:  the  annealing  profile 
at  this  time  differs  very  slightly  from  the  implanted  profile.  Near  the 
profile  peak,  the  boron  has  precipitated  and  it  is  immobile.  The  remain¬ 
ing  boron,  which  is  composed  mainly  of  BV-pairs,  has  diffused  very 
slightly  in  this  short  time.  Later  in  time,  annealing  proceeds  with  two 
mechanisms.  One  mechanism  is  thr  conversion  of  BV-pairs  into  electri¬ 
cally  active  boron  accompanied  by  redistribution  as  described  in  the 
first  canonical  case.  The  other  mechanism  is  the  dissolution  of  precip¬ 
itated  boron.  The  dissolution  takes  place  because  the  redistribution  of 
boron  removes  boron  from  the  areas  of  high  concentration,  thus  upsetting 
the  equilibrium  concentrations  of  precipitated  and  nonprecipitated  boron 
with  respect  to  solid  solubility.  This  nonequilibrium  situation  is,  in 
turn,  corrected  by  the  dissolution  of  precipitates.  Whereas  the  B-BV  pair 
reaction  reaches  equilibrium  conditions  in  the  course  of  annealing,  dis¬ 
solution  a  slow  process  at  900 °C  and,  after  35  min  of  annealing,  there 

1-  still  a  large  fraction  of  boron  precipitates.  This  situation  is  de¬ 
picted  in  the  insert  to  the  left  in  Fig.  12b.  The  boron  precipitate  re¬ 
mains  essentially  unchanged,  however,  the  fraction  of  BV-pairs  diffused 
as  they  converted  into  substitutional  boron.  Hence,  the  total  boron  con¬ 
centration  profile  exhibits  a  shoulder-like  structure  in  agreement  with 
the  experiment.  Furthermore,  we  can  also  predict  that  subsequent  anneal¬ 
ing  will  produce  diffusion  of  the  tail  of  the  profile  at  ordinary  diffu¬ 

sion  rate  because  the  B-BV  pair  reaction  reached  equilibrium  conditions, 
and  the  electrical  activation  will  be  extremely  slow  because  equilibrium 
with  respect  to  solid  solubility  is  only  slightly  perturbed  by  the  slow 
diffusion  at  900°c.  To  elaborate  further  on  this  subject,  we  can  con¬ 
trast  this  case  with  a  higher  annealing  temperature  case.  As  the  anneal¬ 
ing  temperature  is  increased,  the  solid  solubility  limit  increases  and 
diffusion  becomes  more  important.  Consequently,  the  dissolution  of  boron 
precipitates  will  occur  more  readily  and  a  higher  rate  of  electrical  ac¬ 
tivation  results.  For  instance,  at  1100°C,  the  electrical  activity  ap¬ 
proaches  100^  after  35  minutes  of  annealing.  In  summary,  annealing  of 
high  dose  implants  occur  as  a  result  of  two  mechanisms  .  The  boron-BV 
pair  reaction  governing  the  redistribution  behavior  and  the  dissolution 
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of  boron  precipitates  governing  the  annealing  time.  We  model  the  pre¬ 
cipitation  and  dissolution  kinetics  by  assuming  again  a  rate  propor¬ 
tional  to  the  departure  from  equilibrium  as  the  driving  force  and  a 
time  constant.  Thus,  the  precipitation  rate  becomes: 

C  -  Solid  Solubility  Limit  (SSL) 

Precipitation  Rate  =  +  -  (11 

T 

prec 


and  likewise  the  dissolution  rate  results  in  a  more  complicated  rela- 
t  ion  : 


The  upper  branch  accounts  for  the  situation  depicted  in  Fig.  13a,  in 
which  the  fraction  of  precipitates  below  the  solid  solubility  limit 
(SSL)  indicated  by  the  bracket  should  decay  to  zero.  Similarly,  the 


Fig.  13. 


THE  DISSOLUTION  RATE  OF  PRECIPITATES. 


lower  branch  represents  the  case  in  Fig.  13b  when  the  total  boron  con¬ 
centration  is  below  the  solid  solubility  limit.  In  this  case,  if  there 
is  a  precipitate,  it  should  decay  spontaneously  to  zero. 


By  inspection  of  these  equations,  we  can  verify  that  the  rates 
reduce  to  zero  at  equilibrium. 

7 .  Low  Temperature  Annealing 

We  will  examine  now  the  low  temperature  annealing  cases  in 

which  'metastable  boron  precipitates  ’  form  as  a  result  of  the  trapping 

of  boron  by  dislocation  dipoles.  The  canonical  case  representative  of 

14  2 

this  group  is  the  10  ions/cm  dose  implant  annealed  at  800 °c .  In  this 

„  18  ,3 

example,  the  boron  peak  concentration  is  8  x  10  atoms/cm  ,  which  is 
19  3 

below  10  atoms/cm  ,  the  solid  solubility  limit  at  800°C.  However,  the 
experimental  profile  that  results  after  35  minutes  of  annealing  shows 
again  the  shoulder-like  structure  in  Fig.  14a.  The  electrical  carrier 
profile  show's  a  behavior  very  similar  to  the  high  dose  cases  with  the 
presence  of  true  precipitates.  In  essence,  it  appears  that  there  is 
again  a  threshold  level  of  concentration,  resembling  a  solid  solubility 
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Fig.  14.  THE  ELECTRICAL  ACTIVATION  AND  REDISTRIBU¬ 
TION  BEHAVIOR  OF  BORON  IN  LOW  ANNEALING  TEMPERATURE 
CASES  . 


limit,  above  which  concentration  the  boron  becomes  immobile  and  electri¬ 
cally  inactive.  The  boron  in  this  case  is  certainly  not  a  precipitate 
in  the  conventional  sense.  An  analysis  of  transmission  electron  micro¬ 
graphs  and  electrical  activation  results  indicate  that  this  boron  is 
trapped  by  dislocation  dipoles  that  form  and  evolve  during  annealing.  A 
more  detailed  description  of  this  issue  was  given  in  a  previous  section. 
Our  objective  here  is  to  describe  the  role  and  participation  of  the 
trapped  boron  in  the  low  temperature  annealing.  At  the  outset  of  the 
anneal,  the  large  concentration  of  vacancies  generated  by  the  implanta¬ 
tion  coalesce  into  vacancy  clusters  and  loops,  as  indicated  by  the  5 
minute  transmission  electron  micrograph  in  Fig.  9b.  We  assume  that  the 
low  electrical  activity  is  due  to  a  large  concentration  of  BV-pairs,  as 
shown  in  Fig.  14b.  As  annealing  proceeds,  the  micrograph  corresponding 
to  in  minutes  of  annealing  (Fig.  9b)  indicates  that  the  vacancy  clusters 
break  up  and  two  events  take  place.  One  is  the  migration  of  vacancies 
in  the  bulk  and  to  the  surface.  The  other  is  the  formation  of  disloca- 
t  ion  di Doles  in  the  substrate.  The  consequences  of  this  annealing  stage 
is  the  attainment  of  equilibrium  concentration  for  the  vacancies  and  the 
trapping  of  boron  atoms  by  the  dislocation  dipoles.  This  is  depicted  in 
Fig.  14b,  and  we  can  see  that  the  trapped  boron,  immobile  and  electri¬ 
cally  inactive,  is  formed  at  the  expense  of  the  BV-pairs.  Apparently, 
the  formation  of  dislocation  dipoles  is  very  fast  in  time,  and  the  boron 
immobilized  by  these  defects  is  the  fraction  of  boron  that  exceeds  a 
critical  concentration.  From  these  premises,  we  model  the  trapping  of 
boron  as  a  chemical  reaction  with  a  reaction  rate  represented  by: 


CBV  Cdd 
Tdd 


(13) 


■vtiere  C',,  represents  the  aforementioned  critical  concentration, 
dd 

There  are  other  evidences  of  the  existence  of  this  critical 
concentration  level.  For  instance,  dislocation  dipoles  have  only  been 
observed  in  high  concentration  diffusions  and  not  at  lower  concentra¬ 
tions.  Furthermore,  since  dislocation  dipoles  form  when  they  are  about 
li)0  long  (Bicknell  et  al  [17]),  we  can  estimate  the  concentration 
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Such  a  calculation 


level  for  which  two  boron  atoms  are  10O  \  appart . 

18  'A 

yields  a  boron  concentration  of  10  atoms/cm  ,  which  is  consistent  with 
the  observed  threshold  . 

Once  the  boron  has  been  trapped  by  the  dislocation  dipoles , 
further  annealing  proceeds  with  two  distinct  mechanisms.  One  is  the  B- 
BV  pair  reaction,  as  in  the  other  cases  that  we  already  examined.  The 
other  mechanism  is  the  release  of  trapped  boron  atoms  from  the  disloca¬ 
tion  dipoles  as  they  anneal  out.  The  transmission  electron  micrograph 
at  15  minutes  of  annealing  indicates  that  the  dislocation  dipoles  have 
decreased  in  number  and  increased  somewhat  in  length.  Hence,  wo  would 
expect  a  release  of  trapped  boron  atoms.  We  do  not  have  clear  evidence 
of  how  the  trapped  boron  becomes  active.  However,  the  faster  ci  if  fusion, 
inferred  from  Table  2  on  the  behavior  of  the  annealed  distribution,  ap¬ 
pears  to  indicate  that  the  trapped  boron  becomes  a  BV-pair,  the  fast 
diffusart.  It  is  conceivable  that  a  vacancy  migrates  to  the  neighbor¬ 
hood  of  a  boron  atom  imbedded  in  the  strain  field  of  a  dislocation  di¬ 
pole.  '.Mien  the  dislocation  dipole  is  annealing  out,  the  strain  field 
is  weakened  and  the  vacancy  may  approach  the  boron  atom  to  form  a  boron 
vacancy  pal’-.  It  is  also  possible  that  the  dislocation  dipoles  emit 
vacancies  as  they  anneal  out,  hence  increasing  locally  the  probability 
for  the  formation  of  BV-pairs.  The  global  features  characteristic  of 
this  case  are  depicted  in  Fig.  11b.  The  inserts  below  the  graph  repre¬ 
sent  schematically  the  evolution  of  defects,  among  which  are  the  dislo¬ 
cation  dipoles,  capable  of  trapping  boron  atoms  in  their  strain  fields. 
In  the  electrical  activity  graph,  we  show,  using  a  bar  graph  represen¬ 
tation,  the  fractions  of  immobile  boron,  BV-pairs,  and  electrically  ac¬ 
tive  boron .  The  arrows  connecting  these  fractions  represent  the  contin¬ 
uous  release'  of  trapped  boron  from  the  dislocation  dipoles  astheyanneal 
out  and  the  conversion  of  BV-pairs  into  active  boron.  The  inserts  above 
the  electrical  activation  graph  indicate  the  evolution  of  the  redistri¬ 
bution  behavior  that  results  from  the  presence  of  immobile  boron,  BV- 
pairs  (fast  diffusant),  and  substitutional  boron  (slow  diffusant).  The 
shoulder-like  structure  at  annealing  times  of  15  and  55  minutes  are  con¬ 
sequences  of  the  immobility  of  the  trapped  boron  near  the  profile  peal; 
and  the  mobility  of  the  BV-pairs  responsible  for  the  faster  diffusion  "f 


15 


the  profile's  tail.  We  can  appreciate  at  +  his  point  that  the  major  dif¬ 
ference  bet. seen  this  present  case  and  the  high  dose  case  originates  from 
tlie  fact  that,  whereas  the  dissolution  of  boron  precipitates  in  the  high 
dose  case  is  a  consequence  of  the  removal  of  boron  from  the  areas  of  high 
concentration,  in  the  low  annealing  temperature  case  the  release  of  boron 
i'  the  result  of  the  annealing  of  dislocation  dipoles.  The  release  of 
boron  from  the  strain  fields  associated  with  the  dislocation  dipoles  is 
a  complicated  matter.  The  individual  strain  fields  may  overlap  and  pro¬ 
duct!  a  nonlinear  relation  between  the  number  of  boron  atoms  associated 
with  a  dislocation  dipole  and  the  dislocation  dipole  density.  However, 
as  we  will  see,  despite  the  complexity  in  the  real  situation,  the  simple 
exponential  release  rate  in  Eq .  (11)  is  capable  of  capturing  the  key  fea¬ 
tures  under  consideration: 


B  Release  Rate  =  — ti apped  (11) 

Tdd 

This  equation  is  mathematically  identical  to  the  expression  for  the  dis¬ 
solution  of  the  precipitate  when  the  boron  concentration  is  below  the 
solid  solubility  limit.  Therefore,  Eq .  (12)  can  be  used  inthehigh  dose 
as  well  as  in  the  low  annealing  temperature  cases,  provided  the  appro¬ 
priate  parameters  are  substituted. 

Let  us  examine  now  the  global  features  of  the  anneal  when  im¬ 
plantation  dose  or  annealing  temperature  are  increased  with  respect  to 
the  canonical  case  that  we  just  studied.  An  increase  in  dose  produces 
a  proportional  increase  in  the  dislocation  dipole  density.  The  time 
evolution  of  the  defects  is  depicted  in  Fig.  Bb .  The  principal  differ¬ 
ence  with  the  canonical  case  is  the  formation  of  residual  damage  from 
ttie  annealing  of  the  dislocation  dipoles.  In  contrast  with  the  dipoles, 
the  residual  damage  does  not  trap  boron  atoms.  Annealing  of  the  higher 
dose  case  at  the  same  annealing  temperature  of  sno°c  requires  a  longer 
annealing  time  because  a  larger  fraction  of  boron  is  trapped  by  the  de¬ 
fects.  Tlie  other  situation  of  interest  is  when  the  annealing  tempera¬ 
ture  is  increased.  At  a  h igher  temperature,  annealing  becomes  more 
For  instance,  Fig.  H  shows  that  the  dislocation  dipole  density 
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after  5  minutes  of  annealing  at  900°C  is  equivalent  to  35  minutes  of 
annealing  at  800°C.  Presumably,  at  higher  temperature,  the  sequence  of 
events  in  the  800 °c  anneal  depicted  in  the  micrographs  in  Fig.  8b  at  5, 
10,  and  15  minutes  occur  at  earlier  times  in  the  900 °C  anneal.  Conse¬ 
quently,  annealing  in  this  case  is  not  dominated  by  the  time  required 
to  release  the  trapped  boron  but,  instead,  is  governed  by  the  B-BV  pair 
reaction  time  constant.  This  is  indeed  the  situation  idealized  by  the 
first  canonical  case  in  which  the  presence  of  boron  trapped  by  the  dis¬ 
location  dipoles  is  assumed  to  be  negligible. 

8 .  Summary 

We  describe  in  this  chapter  the  annealing  behavior  of  boron 

that  is  ion  implanted  at  room  temperature  into  silicon,  in  the  dose 
,  ,.14  16  ,2 

range  of  10  to  10  ions/cm  and  annealed  in  the  temperature  range  of 
80< i °C  to  1000 °C.  This  wide  range  of  implantation  dose  and  annealing 
temperatures  include  from  the  typical  annealing  cases  of  ion  implanta¬ 
tion  in  the  fabrication  of  devices  to  the  unconventional  situations  in 
which  the  annealing  is  dominated  by  boron  precipitates  (high  dose  im¬ 
plantation)  and/or  boron  trapped  by  dislocation  dipoles  (low  temperature 
a nnea  Ling)  . 

To  model  the  redistribution  and  electrical  activation  of  boron, 
we  propose  a  model  composed  of  substitutional  boron  (electrically  active 
and  slow  diffusant),  boron-vacancy  pair  (electrically  inactive  and  fast 
diffusant),  and  the  positively  charged  vacancy.  The  overall  diffusion 
of  boron  is  then  composed  of  a  slow  and  electrically  active  fraction, 
the  substitutional  boron,  and  a  fast  and  electrically  inactive  fraction, 
the  boron  vacancy  pair.  The  diffusive  and  electrical  attributes  of  boron 
and  the  BV-pair  enable  the  overall  diffusion  of  boron  to  assume  the  di¬ 
chotomous  characteristics  that  pertain  to  the  prevalent  specie  in  the 
diffusion.  For  instance,  the  implanted  boron  react  with  nearby  vacancies 
to  form  BV-pairs.  This  dominant  population  of  BV-pairs  is  responsible 
for  the  low  electrical  activity  and  enhanced  diffusion  at  the  outset  of 
the  anneal.  During  the  anneal,  the  BV-pairs  break  up  spontaneously  into 
its  constituents  until  the  thermodynamic  equilibrium  is  reached.  At  equi¬ 
librium,  the  opposing  rates  in  the  reaction  cancell  and  the  remaining 
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fraction  of  BV-pairs  is  small  compared  to  the  population  of  substitu¬ 
tional  boron.  For  this  reason,  the  diffusion  is  slow  and  the  electri¬ 
cal  activity  is  high  at  the  completion  of  the  anneal.  This  annealing 
behavior  is  characteristic  of  the  typical  annealing  cases.  In  order 
to  include  the  precipitation  effects  of  high  dose  implants,  and  the  low 
annealing  temperature  anomalies  caused  by  the  trapping  of  boron  by  dis¬ 
location  dipoles  ,  it  is  necessary  to  expand  this  basic  model  to  include 
the  immobile  boron  and  additional  reaction-kinetic  terms  .  In  the  next 
chapter,  we  will  utilize  these  premises  to  formulate  a  mathematical 
model  for  the  annealing  of  ion  implanted  boron  at  room  temperature, 
under  the  described  conditions  of  implantation  dose  and  annealing  tem¬ 
pera  ture  . 
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Chapter  II 


THE  ANNEALING  MODEL 


In  this  chapter,  we  will  develop  the  mathematical  model  for  predict¬ 
ing  the  annealing  behavior  of  ion  implanted  boron  into  silicon.  Because 
annealing  is  a  very  complex  problem,  the  first  chapter  was  devoted  to  the 
description  and  definition  of  the  problem,  to  the  analysis  of  relevant 
experiments  pertinent  to  annealing.  From  these  experiments,  we  inferred 
the  key  attributes  of  the  annealing  problem.  The  results  of  these  infer¬ 
ences  are  summarized  in  Table  1.  In  the  left-hand  side  of  the  table,  we 
show  the  species  and  their  electrical  and  diffusive  characteristics.  To 
the  right,  we  indicate  the  reactions  the  species  participate  in,  .-.hich  an¬ 
nealing  behavior  a  particular  reaction  determines,  and  when  this  reaction 
orevail'-.  It  is  apparent  that  the  kinetics  of  the  dominant  reaction  and 
the  electrical  and  diffusive  attributes  of  reactants  and  products  deter¬ 
mine  the  electrical  activation  and  the  redistribution  features,  respec¬ 
tively.  For  thi~  reason,  it  is  important  to  relate  the  reaction  that  pre¬ 
vails  v.  i  t)i  the  correspond  ing  implantation  and  annealing  conditions.  This 
•ross -re  ference  is  provided  by  Table  3,  which  is  repeated  here  from  the 
proved ing  chapter . 
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KKY  ATTRIHl'TES  IN  Till:  ANNEALING  PKOHIJvM 


A. 


The  Mathematical  Model 


The  quantitative  description  of  the  ideas  summarized  in  Table  1  is 
represented  below  by  Eq .  (15).  In  essence  ,  this  set  of  equations  performs 
the  bookkeeping  of  each  of  the  species  in  time  and  space  with  regard  to 
entries  such  as  diffusion,  chemical  reactions,  etc. 
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DIFFUSION  REACTIONS 


iPREC)  (DISS) 
TRAP)  "  (REL  f 


(15a) 


(15b) 


(15c) 


(15d) 


where 

2 

C  =  atomic  concentration  in  #/cm“ 

2 

D  =  diffusion  coefficient  in  cm  /sec 


and  the  subscripts  denote 


B  =  substitutional  boron,  electrically  active 
V+  =  positively  charged  vacancy 
BV  =  boron  vacancy  pair 

Bi  =  immobile  boron,  a  boron  precipitate  in  the  high  dose  implanta¬ 
tion  cases  or  a  boron  trapped  by  dislocation  dipoles  in  the  low 
temperature  annealing  cases 
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The  brackets  in  Kq .  (12)  denote  alternative  choices.  The  upper  selec¬ 
tion  in  the  brackets  corresponds  to  high  dose  implantation  cases,  and 
the  lower  selection  pertains  to  low  annealing  temperatures.  The  ana¬ 
lytical  expression  of  the  various  reactions  in  Eq .  (15)  are  summarized 
in  Eq  .  ( 1(3  )  . 
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k  =  equilibrium  constant,  a  function  of  temperature  only  (cm  ) 
o 

r  -  reaction  time  constant 


■lability  limit.,  function  of  temperature 
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The  reasons  for  representing  boron  precipitates  (high  close  implants) 
and  boron  trapped  by  dislocation  dipoles  (low  annealing  temperatures)  as 
one  -ingle  speries  in  Eq .  (lb)  are  threefold.  First,  we  recall  that  boron 
precipitate's  and  trapped  boron  are  both  electrically  inactive  and  immobile, 
.second,  a  c<>i:.pari -on  of  Kq  .  (lot))  with  (lbc)  and  Kq .  (Kid)  with  (16e)  re¬ 
veals  that,  with  t ho  exception  of  the  parameters,  the  analytical  expres¬ 
sion;  are  respectively  identical  .  Third,  the  only  case  in  which  boron 

l(i 

precipitate  arid  trapped  boron  are  simultaneously  present  is  when  the  W 
•> 

ions /cm  dose  boron  implant  is  annealed  at  XOO °C .  In  this  case,  the  elec¬ 
trical  activation  of  boron  due  to  the  release  of  boron  trapped  by  dislo¬ 
cation  dipoles  is  so  slow  that  the  distinction  between  the  two  forms  of 
boron  becomes  unnecessary.  In  th i -  case  also,  the  diffusion  is  so  slight 


transmission  electron  micrographs  of  samples  annealed  at  low  tempera¬ 
tures  for  short  times  (Figs.  9  and  10).  However,  for  the  case  of  the 
precipitates  present  in  high  dose  implants ,  this  is  strictly  an  assump¬ 
tion.  An  alternative  approach  of  assigning  fractions  of  the  as-implanted 
boron  profiles  to  the  BV-pairs  and  immobile  boron  could  be  as  adequate. 
Transmission  electron  micrographs  also  indicate  that  the  high  vacancy 
concentration  in  Fig.  11  will  condense  first  into  vacancy  clusters  and 
vacancy  dislocation  loops  imbedded  in  a  background  of  vacancies  at  equi¬ 
librium  concentration.  Then,  these  defects  will  anneal,  transform,  and 
release  vacancies  very  rapidly.  For  this  reason,  we  use  the  equilibrium 
concentration  of  positive  vacancies  as  the  initial  condition.  In  most 
instances,  the  release  of  vacancies  from  the  annealing  defects  has  1  ittle 
overlap  in  time  with  the  boron-BV-pa ir  reaction.  Therefore,  the  effects 
re  negligible.  And,  when  the  vacancies  released  from  the  vacancy  clusters 
and  dislocation  loops  do  interact  with  the  annealing,  we  use  the  genera¬ 
tion  term  in  the  vacancy  equation  to  model  a  pulse  of  vacancies  that  is 
distributed  in  space  in  accordance  with  the  energy  deposition  profile. 

In  summary,  the  initial  conditions  are: 


We  recall  that  the  primary  input  to  the  model  is  the  history  of  the  sub¬ 
strate  summarized  by  two  distributions:  the  as-implanted  boron  profile 
and  the  distribution  of  energy  deposited  into  nuclear  processes  by  the 
boron  ions.  We  can  appreciate  now  that,  in  this  implementation  of  the 
annealing  problem,  the  as-implanted  B  profile  becomes  an  initial  condi¬ 
tion,  while  the  energy  deposition  profile  becomes  the  space  dependence 
of  the  vacancy  generation  term. 


(a)  a  set  of  parameters  as  estimated  in  the  preceding  section 

(b)  a  set  of  initial  conditions,  one  for  each  of  the  indepen¬ 
dent  variables: 


Naturally,  the  initial  conditions  are  particular  to  a  specific  problem. 
In  this  section,  our  objective  is  to  elaborate  on  those  initial  condi¬ 
tions  pertinent  to  the  annealing  of  B  in  Si. 

The  detailed  discussion  with  regard  to  the  overall  annealing  be¬ 
havior  can  be  found  in  Chapter  I.  In  that  discussion,  special  attention 
is  devoted  to  the  state  of  the  boron-silicon  system  at  the  outset  of  the 
anneal.  In  particular,  the  initial  state  of  the  silicon  substrate  con¬ 
taining  the  as-implanted  boron  atoms  is  carefully  analyzed.  For  this 

reason,  we  can  formulate  the  initial  conditions  of  interest  by  simply 

recalling  the  results  in  Chapter  I. 

The  as-implanted  profile  as  determined  by  secondary  ion  mass  spec¬ 
trometry  is  shown  in  Fig.  11.  The  vacancy  concentration  profile  is  cal¬ 
culated  from  the  Frenkel -pairs  produced  by  the  70  KeV  incident  boron 
ions  as  they  deposit  energy  into  atomic  processes  [38]  along  their  tra¬ 
jectories.  In  what  follows,  we  will  make  the  conservative  assumption 

that  as  the  boron  atoms  come  to  rest  they  will  all  combine  with  nearby 

vacancies  to  form  BV-pairs,  which  are  electrically  inactive.  Hence,  at 
the  outset  of  the  anneal,  the  electrical  activity  is  zero.  Although  a 
nonzero  activity  is  more  reasonable,  it  will  also  require  an  initial 
profile  for  the  electrically  active  boron.  At  present,  we  do  not  w'ish 
to  introduce  this  unknown.  In  other  words,  the  initial  profiles  of 
boron  and  BV-pair  are  the  null  and  the  as-implanted  profiles,  respec¬ 
tively.  With  regard  to  the  immobile  boron  present  in  high  dose  implant 
and  low  annealing  temperature  cases,  we  assume  that  the  initial  concen¬ 
tration  profile  is  identically  zero.  This  is  in  agreement  with  the 


(hi 


For  these  reasons,  the  only  requirement  for  the  value  of  this  time  con¬ 
stant  is  for  it  to  be  small.  We  choose  both  parameters  to  be  1  sec  ar¬ 
bitrarily  . 

10.  Time  Constant  Associated  with  the  Dissolution  of  Precipitates 
and  the  Release  of  Boron  from  Dislocation  Dipoles 

In  the  higli  dose  implant  cases,  the  dissolution  is  so  slow 

8 

compared  to  the  35  minute  anneal  that  xnTQC  is  chosen  to  be  10  sec, 

Dloo 

an  arbitrarily  large  value.  In  the  low  annealing  temperature  cases,  we 

estimate  x  from  the  electrical  activation  data  by  Seidel  and  MacRae 

REL 

[18].  According  to  these  authors,  the  time  to  reach  90^  electrical  ac- 
„  3  5 

tivity  at  800  c  are  B  x  10  and  1.8  10  sec  for  implantation  doses  of 

1011  and  1015  ions  /cm"" ,  respectively.  We  use  the  values  of  2.4  10  and 
5 

3.6  10  sec  in  the  calculation. 

In  summary,  Table  6  presents  the  global  features  of  the  an¬ 
nealing  problem  by  displaying  the  relation  of  implantation  dose  and  an¬ 
nealing  temperatures  to  the  parameters  in  the  model,  hence  to  the  vari¬ 
ous  mechanisms  that  prevail  during  annealing.  The  parameters  associated 
with  the  annealing  of  damage,  high  dose  implants,  and  low  annealing  tem¬ 
peratures  are  estimated  from  a  more  limited  amount  of  data.  Functionally 
they  account  for  the  formation  of  the  immobile  and  electrically  inactive 
phase  of  boron  and  the  subsequent  transformation  into  BV-pairs.  Here¬ 
after,  the  behavior  of  the  BV-pairs  is  governed  by  the  set  of  parameters 
in  Tabic  5.  This  other  subset  of  parameters  describe  the  redistribution 
and  electrical  activation  of  ion  implanted  boron.  They  describe  ordinary 
diffusion  with  high  dose  anomalies.  They  describe  the  enhanced  diffusion 
of  boron  and  subsequent  reduction  to  ordinary  diffusion.  And,  they  de¬ 
scribe  the  enhancement  in  diffusion  produced  by  proton  implants.  This 
subset  of  parameters  in  the  model  encompass  the  main  issues  in  the  an¬ 
nealing  problem  and  their  importance  cannot  be  overstated. 

C .  The  Initial  Conditions 

To  solve  the  initial  value  problem  described  by  the  set  of  coupled 
partial  differential  equations,  we  need  the  following: 


8 .  Threshold  Concentration  I;evel  for  the  Formation  of  Dislocation 
Dipoles 

We  use  this  parameter  to  represent  the  trapping  of  boron  by 
dislocation  dipoles.  It  is  obvious  that,  functionally,  this  threshold 
level  will  produce  the  desired  fraction  of  electrically  inactive  and  im¬ 
mobile  boron  that  we  propose  to  be  trapped  by  dislocation  dipoles.  Phys¬ 
ically,  this  parameter  represents  the  composite  effects  of  a  sequence  of 
complex  events  depicted  by  the  transmission  electron  micrographs;  viz., 
the  formation  of  dislocation  dipoles  and  the  trapping  of  boron.  As  a 

consequence,  C  cannot  be  obtained  from  simple  physical  arguments. 
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The  value  of  1  X  10  cm  for  C  at  800 °C  is  first  obtained  from  the 
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inspection  of  Fig.  4  (10  dose)  and  then  varied  to  optimize  the  fit  to 
experimental  data.  At  higher  temperatures,  in  the  lower  left-hand  side 
of  Table  6d ,  we  can  infer  from  the  transmission  electron  micrographs  at 
900 °C  that  the  formation  of  dislocation  dipoles  and  their  subsequent  an¬ 
nealing  occurs  so  rapidly  that  the  effects  on  annealing  are  not  observ¬ 
able.  It  is  for  this  reason  that  the  corresponding  values  of  C  are 
not  estimated.  At  the  right-hand  side  of  Table  Gd ,  we  assume  that  all 
the  immobile  boron  are  precipitates  that  arise  from  exceeding  the  SSL. 
For  this  reason,  we  use  the  values  of  the  solid  solubility  limit  for  C 

in  this  region  of  the  table.  Lastly,  TEM  indicates  the  presence  of  dis- 

1 5 

location  dipoles  during  the  900°c  annealing  in  the  case  of  the  10  ions/ 
2 

cm"  dose  implant.  And,  we  find  that,  by  incorporating  the  effects  of 

dislocation  dipoles  in  the  conventional  annealing  model,  we  can  improve 

slightly  the  fit  between  the  calculated  and  the  experimental  profiles. 

.  19  -3  .  _  . 

We  use  the  value  of  o  '  10  cm  for  C  in  this  case. 

9 .  Time  Constants  Associated  with  Precipitation  and  the  Trapping 
of  Boron  by  Dislocation  Dipoles 

In  the  low  annealing  temperature  cases,  the  sequence  of  trans¬ 
mission  electron  micrographs  show  that  the  trapping  of  boron  by  disloca¬ 
tion  dipoles  occurs  early  during  annealing.  In  the  high  dose  implant 
cases,  we  assume  that  the  precipitation  of  boron  takes  place  rapidly.* 

* 

It  is  possible  that  precipitation  has  occurred  prior  to  annealing. 


In  summary,  with  the  exception  of  t  ,  which  is  assumed  to 

V 

be  constant,  all  the  other  parameters  are  functions  of  temperature  only 
To  emphasize  this  dependence,  we  use  the  following  representation.  The 
dotted  lines  within  the  matrix  (a)  in  Table  6  enclose  sets  of  values 
that  are  invariant  with  implantation  ciose,  but  dependent  on  the  anneal¬ 
ing  temperature.  As  mentioned  previously,  the  parameters  that  follow 
hereafter  are  associated  with  anomalies  arising  from  the  annealing  of 
implantation  damage,  high  dose  implants,  and  low  annealing  temperatures 


6.  Damage  Annealing  Time  Constant 

This  parameter  characterizes  the  time  evolution  of  vacancy 
clusters  during  annealing.  Figures  9  and  10  show  the  presence  of  these 
vacancy  sources  in  the  transmission  electron  micrographs  corresponding 
tc  3  minutes  of  annealing  at  800°C.  At  800°C  and  900°C,  because  of  the 
longer  lifetime  of  BV-pairs  10  min),  the  release  of  vacancies  from 
vacancy  clusters  (  :  10  min)  does  not  alter  the  boron-BV-pa ir  reaction 
significantly.  On  the  other  hand,  at  1000°C  the  lifetime  of  BV-pairs 
is  only  110  sec.  Consequently,  the  simultaneous  generation  of  vacancies 
slows  down  the  conversion  of  BV-pairs  into  active  boron,  prolonging  lo¬ 
cally  the  enhanced  diffusion.  This  effect  is  only  observable  at  high 
temperatures  and  short  annealing  times  because,  at  longer  annealing 
times,  it  is  masked  by  ordinary  diffusion.  In  our  calculations,  we  use 

the  value  of  40  sec  for  t  at  1000°C  (see  Table  6b). 

DAM 

At  lower  temperatures,  we  can  deduce  from  transmiss ion  electron 

micrographs  and  the  lifetime  of  BV-pairs  that  the  bounds  of  t  are 

DAM 

larger  than  40  sec  and  smaller  than  550  sec,  respectively. 


7 .  Solid  Solubility  Limit  (SSL) 

The  solid  solubility  limit  is  a  function  of  temperature  only. 
Table  6c  shows  the  values  we  use  in  the  calculation.  These  values  are 
comparable  to  the  experimental  values  by  Vick  et  al  [431. 
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Table  8 


THE  ESTIMATION  OF  THE  LIFETIME  OF  BV-PAIRS 


x  Measured 

from  the 
Annealing  of 
10^4  Dose 
Implants 
(pm) 

T  Estimated 
from  Dbv 
Used 
(sec) 

T  Estimated 
from 

Enhanced 

Diffusion 
Coefficient 
(sec ) 

T 

Used 

(sec) 

800  °C 

0.110 

3000 

— 

2780 

900  °C 

0.079 

230 

1150 

550 

1000 °c 

0.180 

270 

540 

140 

repeat  the  calculation  using  the  enhanced  diffusion  coefficient  deter¬ 
mined  by  Hofker  et  al ;  the  results  of  this  second  estimation  are  shown 
in  the  third  column.  An  alternative  way  to  obtain  the  value  of  t  is 
to  solve  the  set  of  diffusion  equations  with  initial  conditions  for  an 
ordinary  diffusion.  Then,  the  values  of  r  are  chosen  to  fit  experi¬ 
mental  ordinary  diffusion  profiles.  The  values  of  t  we  use  in  the 
calculations  are  presented  in  the  fourth  column  of  Table  9. 

5 .  The  Lifetime  of  Vacancies 

The  lifetime  of  vacancies  is  estimated  from  the  diffusion 
length  of  vacancies.  We  use  the  diffusion  length  of  1.78  X  10  cm  at 
750 °C ,  after  Tsuchimoto  et  al  ri2l.  The  corresponding  value  of  com¬ 

puted  from  this  diffusion  length  and  the  diffusion  coefficient  estimated 

-  5  -4 

previously  is  8.3  10  sec,  approximately  10  sec.  It  is  apparent  that 
r  is  much  smaller  than  r.  Consequently,  provided  that  this  strong 
inequality  is  preserved,  the  annealing  calculations  will  be  rather  in¬ 
sensitive  to  the  actual  value  of  r  .  For  this  reason,  the  lack  of  more 
data  on  the  diffusion  length  of  vacancies  does  not  create  a  problem,  and 
we  can  use  the  value  of  10  *  sec  for  r^,  at  all  temperatures .  It  is 
evident  that,  if  more  data  were  available,  then  we  could  have  estimated 
a  vacancy  lifetime  that  is  temperature  dependent. 


03 


I 


( 


where  Ey+  =  0*35  eV  above  the  valence  band  edge.  The  Fermi  level  is 
calculated  using  approximations  to  the  Fermi-Dirac  integral  of  order  1/2 
(Blakemore  [41]).  Hence,  doping  degeneracies  are  accounted  for.  It  is 
worthwhile  to  emphasize  that  the  apparent  concentration-dependent  diffu¬ 
sion  coefficient  of  boron  is  achieved  in  this  model  through  the  Fermi 
level.  A  high  substitutional  boron  concentration  shifts  the  Fermi  level 
closer  to  the  valence  band  edge,  thus  increasing  the  concentration  of 
positive  vacancies  (E^  =  0.35  eV  above  the  valence  band).  This,  in 
turn,  increases  the  ratio  of  BV-pairs  to  substitutional  boron,  increas¬ 
ing  the  overall  diffusion  of  boron. 

4 .  t,  The  Lifetime  of  BV-Pairs 

This  parameter  is  estimated  from  the  diffusion  length  associ¬ 
ated  with  the  BV-pairs,  the  fast  diffusing  specie  of  finite  lifetime. 

If  we  assume  that  the  distance  x  travelled  by  the  profile  tail  is 
twice  the  diffusion  length,  then  we  can  solve  for  t  in  the  following 
expression : 


x 


2  ;d 

V  BV 


T 


where  the  distance  x  is  measured  from  the  isochronal  annealing  pro¬ 
files  in  Fig.  5a,  after  Hofker  et  al.  The  values  of  x  at  annealing 
temperatures  in  the  range  of  800°C  to  1000°C  are  shown  in  the  first 
column  of  Table  8.  It  is  appropriate  to  comment  that  the  values  of  x 
thus  obtained  are  more  representative  of  the  enhanced  d if fus ion  at  800°C 
than  at  1000°C.  The  reason  for  this  is  that,  for  a  fixed  annealing  time 
of  35  minutes,  the  contribution  of  ordinary  diffusion  is  less  signifi¬ 
cant  at  lower  temperatures  than  at  higher  temperatures.  Consequently, 
the  estimate  of  x  at  1000°c  is  probably  an  upper  limit.  The  values  of 
x  calculated  in  this  manner  using  the  values  of  D  estimated  previ¬ 
ously  are  shown  in  the  second  column  of  Table  8.  As  a  comparison,  we 
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The  Equilibrium  Constant  in  the  Thermodynamic  Reaction 


2 

At  thermal  equilibrium,  the  concentrations  of  reactants  and 
products  of  the  boron-BV-pa ir  reaction  are  given  by  the  following  ther¬ 
modynamic  equation: 


C 


BV 


k  C  C 
o  B  V 


(4) 


The  equilibrium  constant  k  of  this  reaction  can  be  calculated  from 

o 

Eq .  (4)  as  follow's.  Under  equilibrium  conditions,  such  as  in  an  ordi¬ 
nary  diffusion,  the  electrical  activity  is  near  100’').  Assuming  again 
an  electrical  activity  of  98^  and  rearranging  Eq .  (4),  we  obtain  the 

following  expression  for  k  : 

o 


k 

o 


<CDV/CB> 


equil .  , 


intrinsic , 


(20) 


where  C  +  is  the  equilibrium  concentration  of  positive  vacancies  under 
intrinsic  conditions.  The  value  of  Cy  +  is  obtained  from  Eq .  (8)  with 
the  Fermi  energy  level  at  the  middle  of  the  energy  band  gap.  Conse¬ 
quently,  Cy+  is  a  function  of  temperature  only. 

This  choice  of  Cv+  and  the  method  for  selecting  insure 

that,  when  the  electrical  activity  reaches  98^,  the  overall  diffusion 
coefficient  will  approach  to  the  value  obtained  from  ordinary  diffusion 
experiments  under  intrinsic  conditions.  From  a  different  viewpoint, 
this  choice  of  Cv+  is  responsible  for  the  modelling  of  high  dose  and 
substrate  doping  effects  by  allowing  the  actu  1  concentration  of  posi¬ 
tive  vacancies  to  alter  the  position  of  active  to  inactive  boron  [Ecj . 
(6)1  . 


3 .  Equilibrium  Concentration  of  Positively  Charged  Vacancies 

The  equilibrium  concentration  of  positive  vacancies  C^,+  is 

eq 

dependent  on  the  energy  level  11  ,  ,  the  Fermi  energy  level  E  ;  the 

v  F 

temperature  T  and  the  concentration  of  neutral  vacancies  C  Q ,  which 
is  a  function  of  temperature  only  f36l. 


(it) 


Table  7 


THE  ESTIMATION  OF  D0,r 

BV 


Enhanced 

Diffusion 
Coefficient 
by  Hofker  et  al 
(cm^/sec) 

— 

Dgy  Estimated 
from  the  Dif¬ 
fusion  of 

Annea led 
Profiles 

O 

(cm  /sec) 

Dgy  Estimated 
from  the 
Ordinary 
Diffusion 

Coefficient 

(cm^/sec) 

DBV  Usod 
(cm^/sec ) 

800  °C 

7.5  10-14 

2.51  10-14 

io-14 

900  °C 

1.4  10-14 

1.7  10-14 

7  10"14 

6.5  10“14 

7  X  10-14 

1000 °C 

1.5  10-13 

9  10-13 

8.5  10-13 

-13 

3.56  10 

alternative,  Anderson  and  Gibbons  estimate  the  value  of  D  to  be  5.7  x 
-14  2 

10  cm  /sec  at  750°C  from  proton  enhanced  diffusion  experiments. 

Another  alternative  method  of  estimating  the  value  of  D  is  based  on 

BV 

the  following  relation  between  diffusion  and  electrical  activity: 

D  =  D  (1  -  electrical  activity)  (19) 

exp  BV 

where  D  is  the  diffusion  coefficient  measured  in  an  experiment, 

exp 

hence  the  subindex. 

The  discussion  and  derivation  of  this  relation  is  carried  out 

in  Appendix  B.  To  estimate  D  ,  we  solve  Eq .  (19)  under  equilibrium 

BV 

conditions.  For  this  reason,  we  assume  that  the  electrical  activity  is 

98v  and  we  use  for  D  the  experimental  diffusion  coefficients  mea- 

exp 

sured  under  equilibrium  conditions  by  Hofker  et  al  f2]  and  Kurtz  (see 
Table  1).  The  values  of  thus  estimated  are  shown  in  the  second 

and  third  columns  of  Table  7,  respectively.  In  the  fourth  column,  we 
show  the  values  used  in  the  present  work. 


Table  6 

THE  RELATION  OF  THE  PARAMETERS  TO  THE 
IMPLANTATION  DOSE  AND  ANNEALING  TEMPERATURE 


DB’VDBV’ko'T’TV 


40  sec 


40  sec  40  sec  40  sec 


ions/ cm 

(b) 

(a) 

SSL  (s/cm3) 

0  (s/cm3) 

'  7  X  1018  7  V  1018  7xl018' 


'  4  xlO18  4  v  ID18  4  xlO18' 


*1.1 1020  1.1  1020  1.1  1020' 

I  I 


'3.3  1020  3.3  1020  3.3  1020' 


5  1()19  1.1  X  1020 

(SSL) 

3.3  X 1020 
(SSL) 


tprec  “  tdd  (sec 


1  sec 


2.4  104  3.6  105  i  108* 


Arbitrarily  small 


Arbitrarily  large 


contain  the  preexponential  factors  and  the  activation  energies.  This  set 
of  parameters  is  the  backbone  of  the  annealing  model.  These  parameters 
control  the  redistribution  and  electrical  activation  of  boron.  They  are 
the  dominant  parameters  in  the  ordinary  diffusion  of  boron  and  in  the 
typical  annealing  cases  of  ion  implanted  boron  into  silicon.  In  high 
dose  implantation  and  low  annealing  temperature  cases,  these  parameters 
control  the  redistribution  and  electrical  activation  of  the  boron  aris¬ 
ing  from  the  dissolution  of  boron  precipitates  or  the  boron  released 
from  the  annealing  of  dislocation  dipoles.  As  discussed  previously,  the 
modelling  of  precipitation  and  dissolution  as  well  as  the  trapping  of 
boron  by  dislocation  dipoles  and  the  subsequent  release  require  addi¬ 
tional  parameters  such  as  solid  solubility  limits,  time  constants,  and 
so  on.  These  parameters  are  generally  obtained  from  more  empirical  con¬ 
siderations.  A  summary  of  all  the  parameters  and  their  relation  to  im¬ 
plantation  dose  and  annealing  temperature  is  given  in  Table  6.  In  this 
manner,  Table  6  displays  mainly  the  parameters  associated  with  anomal¬ 
ies,  while  Table  5  contains  the  parameters  controlling  the  main  diffu¬ 
sion  mechanism.  The  data  in  Table  6  is  arranged  in  six  separate  groups. 
Each  group  is  a  matrix  displaying  the  behavior  of  the  parameter  as  im¬ 
plantation  dose  (horizontal)  and  annealing  temperature  (vertical)  are 
varied  . 

We  shall  now  discuss  the  methods  used  in  the  estimation  of  the  pa¬ 
rameters  in  Tables  5  and  6 . 

1 .  Diffusion  Coef f icients--Only  Functions  of  Temperature 

We  assume  that  substitutional  boron  diffuses  by  the  same 
mechanism  as  the  self  diffusion  of  silicon.  Therefore,  we  use  the  self 
diffusion  coefficient  of  silicon  [40]  for  D  .  For  positive  vacancies, 
we  use  the  diffusion  coefficient  proposed  by  Seidel  and  MacRae  [37]. 
Finally,  the  diffusion  coefficient  for  BV-pairs  is  estimated  from  vari¬ 
ous  diffusion  experiments.  Hofker  and  coworkers  measure  the  enhanced 
diffusion  of  ion  implanted  boron  after  35  minutes  of  annealing.  The 
enhanced  diffusion  coefficient  they  obtain  is  shown  in  the  first  column 
in  Table  7.  Since  these  are  average  values,  they  provide  only  lower 

limits  for  the  value  of  D  at  the  respective  temperatures.  As  an 

BV 
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on  a  digital  computer.  The  results  of  the  integration  are  concentra¬ 
tions  of  each  of  the  species  as  functions  of  time  and  space.  Details 
concerning  the  solution  are  given  in  Appendix  A.  In  order  to  evaluate 
the  effectiveness  of  the  model,  we  have  to  compare  the  calculated  re¬ 
sults  with  experimental  data.  In  this  respect,  while  the  electrical 
carrier  concentration  profile  can  be  compared  directly  with  the  calcu¬ 
lated  active  boron  profile  at  corresponding  times,  the  comparison  with 
other  experimental  data  requires  some  algebraic  manipulations.  For  in¬ 
stance,  the  total  boron  concentration  determined  experimentally  using 
Secondary  Ion  Mass  Spectrometry  (SIMS)  is  to  be  compared  with  the  sum¬ 
mation  in  Eq .  (17). 

C  (x,t)  =  C  (x,t)  +  C  (x,t)  +  C  (x,t)  (17) 

15  ,  ,  B  BV  B 1 

tota  1 

And,  the  electrical  activation  curve  is  to  be  compared  with  the  spatial 
integration  of  the  calculated  active  boron  profile  in  Eq .  (18). 

/•+  00 

ea(t)  =  I  C  (x,t)  dx 

•A)  B 

I) .  The  Selection  of  Parameters 

The  mathematical  model  presented  in  the  preceding  section  is  the 
result  of  inferences  drawn  from  the  discussions  in  Chapter  I.  Basically, 
the  complex  annealing  problem  is  broken  into  three  simpler  cases.  Then, 
the  reactions  dominant  in  each  of  the  aforementioned  cases  are  identi¬ 
fied.  The  kinetics  of  these  reactions  are  then  modelled  by  first  order 
approximations  involving  parameters  with  physical  interpretations.  Our 
objective  in  this  section  is  to  estimate  the  numerical  values  of  these 
parameters  . 

In  typical  annealing  cases,  the  problem  is  less  complicated  and  the 
experimental  data  is  more  abundant.  Consequently,  it  is  possible  to  ap¬ 
ply  thermodynamic  and  physical  considerations  in  the  estimation  of  the 
parameters.  Table  5  gives  a  summitry  of  the  estimated  values  of  the  pa¬ 
rameters,  which  are  functions  of  temperature  only.  The  last  two  columns 


that  the  dissolution  of  boron  precipitates  is  negligible  (see  Fig.  )  . 

In  other  cases,  either  boron  precipitate  or  trapped  boron  is  present. 
Consequently,  by  defining  the  immobile  boron  to  represent  either  boron 
precipitate  or  trapped  boron,  we  do  not  lose  generality  and  we  do  gain 
simplicity  in  the  mathematical  formulation  of  the  model. 

In  summary,  we  propose  a  mathematical  model  for  predicting  the  re¬ 
distribution  and  electrical  activation  of  ion  implanted  boron  that  is 
subsequently  annealed.  The  general  model  in  Eq  .  (15)  is  aimed  toward 
solving  the  annealing  problem  under  the  implantation  and  annealing  con¬ 
ditions  in  Table  3.  These  conditions  are  by  no  means  too  restrictive; 
in  fact,  most  applications  fall  into  the  subcategory  of  "typical  anneal¬ 
ing  cases"  in  Table  3.  To  solve  the  simple  annealing  cases,  as  dis¬ 
cussed  in  Chapter  I,  only  three  species  are  required.  Consequently,  Eq . 
(15)  can  be  simplified  by  dropping  the  last  equation  and  the  associated 
coupling  terms  from  the  remaining  equations.  The  importance  of  this  sim¬ 
pler  case  is  threefold.  First,  the  annealing  mechanism  in  this  simpler 
case  is  also  present  in  more  complicated  cases.  Second,  in  these  cases, 
after  35  minutes  of  isothermal  annealing,  the  enhanced  diffusion  reduces 
to  ordinary  diffusion.  Hence,  this  set  of  equations  should  also  apply 
to  ordinary  diffusion  if  appropriate  initial  conditions  are  used.  Third, 
the  enhanced  diffusion  mechanism  at  the  outset  of  the  anneal  is  also  re¬ 
sponsible  for  the  enhancement  in  diffusion  produced  by  proton  implants. 
Consequently,  the  three  species  model  should  be  capable  of  predicting 
ordinary  diffusion,  proton  enhanced  diffusion,  and  the  simpler  annealing 
cases  in  Table  3,  And,  if  the  three  species  model  is  extended  to  include 
boron  precipitate  and  boron  trapped  by  dislocation  dipoles,  then  the  more 
general  model  in  Eq .  (13)  can  also  predict  the  high  dose  implant  anci  low 
annealing  temperature  cases  in  Table  3. 

Let  us  outline  now  the  necessary  steps  for  solving  the  equations  of 
the  annealing  model.  First  of  all,  a  set  of  parameters  is  to  be  estimated 
based  on  thermodynamic  considerations  and  experimental  data.  Then  we  will 
recall  the  discussions  in  Chapter  I  and  establish  the  initial  conditions 
for  the  annealing  problem.  With  these  elements,  we  can  now  integrate  the 
set  of  coupled  partial  differential  equations.  This  is  an  initial  value 
problem,  and  the  solution  is  implemented  using  numeri  •!  analysis  methods 


Chapter  III 


COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  RESULTS 

With  the  choice  of  parameters  and  initial  conditions  described  in 
the  preceding  chapter,  we  can  proceed  to  solve  the  annealing  behavior 
of  ion  implanted  boron  into  silicon.  The  results  of  the  calculations 
corresponding  to  the  cases  in  Table  3  can  be  structured  in  three  canon¬ 
ical  groups: 

(1)  typical  annealing  cases 
(3)  high  dose  cases 

(3)  low  annealing  temperature  cases 

Each  group  has  characteristic  features  that  are  manifest  in  the  equa¬ 
tions,  reaction  kinetic  terms,  and  parameters.  And,  these  characteris¬ 
tics  features  are  consequences  of  special  conditions  present  at  certain 
implantation  dose  and/or  annealing  temperatures.  For  these  reasons,  we 
present  the  results  of  the  calculations  following  the  same  organization 

A .  Typical  Annealing  Cases 

The  result  of  the  calculation  with  the  initial  condition  corre- 

11  ,2 

spending  to  the  10  ions/cm  dose  implant  and  the  parameters  for  the 
900 °c  anneal  is  shown  in  Fig.  15.  The  atomic  concentration  is  repre¬ 
sented  on  the  vertical  logarithmic  scale,  and  the  depth  in  microns  into 
the  silicon  substrate  is  represented  linearly  on  the  horizontal  axis. 
The  as-implanted  profile  determined  experimentally  by  Hofker  et  al  T2l 
is  represented  by  the  dotted  line.  Upon  annealing,  the  BV-pairs  in  this 
profile  diffuse  and  convert  into  substitutional  boron.  Hence ,  a  .sequence 
of  electrically  active  boron  profiles  will  develop  with  time.  Figure  15 
shows  calculated  profiles  after  1,  3,  10,  and  35  minutes  of  annealing 
with  corresponding  electrical  activities  of  10,  28,  65,  and  94  percent. 
The  total  boron  concentration  at  35  minutes  is  represented  by  triangles  ; 
the  fit  to  the  solid  line  representing  the  experimental  SIMS  profile  is 
excellent.  Since  the  ordinary  diffusion  of  boron  at  900 °C  would  only 
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Fig.  15.  THE  ANNEALING  OF  A  10  IONS/CM  DOSE 
B  IMPLANT  AT  900 °C .  The  comparison  of  calcu¬ 
lated  and  experimental  results  . 


modify  the  as-implanted  profile  very  slightly,  the  fit  between  the  cal¬ 
culated  results  and  the  deeper  experimental  profile  is  indicative  of  the 
existence  of  enhanced  diffusion  which  progressively  diminishes  to  the 
ordinary  diffusion  rate  as  the  electrical  activity  approaches  100',.  The 
calculated  results  are  thus  in  good  agreement  with  the  experimental  ob¬ 
servations  of  Hof ker  and  coworkers  [2l. 

15 

Figure  16  shows  the  results  of  the  calculation  for  the  10  ions/ 
cm  dose,  which  produces  an  order  of  magnitude  increase  in  the  impurity 
and  damage  profiles.  We  show  again  the  evolution  of  the  calculated 


B  IMPLANT  AT  900 °C.  The  comparison  of  calcu¬ 


lated  and  experimental  results. 


active  boron  concentration  with  time,  and  we  compare  the  calculated  to 
tal  boron  concentration  represented  by  triangles  with  the  experimental 
SIMS  profile  represented  by  the  solid  line.  Again,  the  fit  is  excellent. 
We  wish  to  comment  that  this  basic  result  can  be  obtained  from  two  al¬ 
ternative  calculations,  by  including  or  excluding  the  trapping  effects 
of  boron  by  dislocation  dipoles,  although  there  is  a  slight  loss  of  ac¬ 
curacy  if  the  effects  of  the  dislocation  dipoles  are  neglected.  Figure 
16  depicts  the  results  of  the  calculation  that  includes  the  trapping  of 
boron  by  defects. 


T 


B  IMPLANT  AT  1000°C.  The  comparison  of  calcu¬ 


lated  and  experimental  results  . 


diffusion.  It  is  for  this  reason  that,  for  short  annealing  times ,  the 
inclusion  of  the  generation  term  in  the  vacancy  equation  contributes  to 
the  improvement  of  the  fit  to  the  experimental  data.  Naturally,  at 
longer  annealing  times,  the  enhanced  diffusion  is  overwhelmed  by  the 
greater  ordinary  diffusion  at  higher  temperatures  and  the  aforementioned 
improvement  becomes  less  noticeable.  We  can  state  from  a  different 
viewpoint  that:  for  the  case  of  a  long  annealing  at  high  temperature, 
the  neglect  of  the  vacancy  generation  term  in  the  calculation  will  not 
cause  significant  errors. 
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The  comparison  of  the  experimental  electrical  activation  of  boron 
versus  time  by  Seidel  and  MacRae  [18]  with  the  calculated  results  from 
some  of  the  previous  examples  is  shown  in  Fig.  19.  The  agreement  is 
good  for  long  times.  For  short  times,  we  attribute  the  difference  to 
our  initial  condition  of  zero  electrical  activity  at  the  outset  of  the 
anneal.  As  discussed  previously,  we  expect  improvements  with  a  nonzero 
initial  condition  for  the  electrical  activity.  The  nonzero  initial 
electrical  activity  is  generally  attributed  to  room-temperature  anneal¬ 
ing  of  the  implanted  silicon  substrate.  We  can  apply  this  idea  to  the 
interpretation  of  the  calculated  electrical  activity  result  in  Fig.  19 
and  to  obtain  an  alternative  set  of  initial  conditions.  For  instance, 
if  t  is  the  instant  in  time  at  which  the  desired  electrical  activity 
is  attained,  then  the  annealing  prior  to  t  is  due  to  room-temperature 
annealing  and  only  the  electrical  activation  of  boron  after  t  is  at¬ 
tributable  to  the  annealing  at  the  elevated  temperature.  In  other  words, 


Fig.  19.  THE  COMPARISON  OF  CALCULATED  ELECTRICAL  AC¬ 
TIVITY  VERSUS  TIME  WITH  EXPERIMENTAL  ISOTHERMAL  RE¬ 
SULTS  BY  SEIDEL  AND  MacRAE.  Dose:  1014  ions/cm2, 
900  . 
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the  experimental  result  should  be  compared  with  the  calculated  result 
shifted  in  time  by  t^  to  correct  for  room  temperature  annealing  ef¬ 
fects.  Or,  equivalently,  a  shift  in  time  would  not  be  required  if  the 

boron  and  BV-pair  profiles  at  t  were  the  initial  conditions  used  in 

o 

the  calculation.  To  explore  the  validity  of  these  arguments,  we  use 

the  pertinent  profiles  at  t  =  2.5  minutes  as  initial  conditions  in 

o 

the  subsequent  calculation.  We  find  that  the  calculated  data  points 
represented  by  'x's,  thus  obtained,  do  approximate  the  experimental 
electrical  activation  curve.  We  also  verify  that  shifting  the  time 
axis  of  the  original  calculation  by  t  produces  equivalent  results. 
Namely,  subtraction  of  2.5  minutes  from  the  abscissa  of  each  original 
data  point  (calculated  using  the  zero  initial  electrical  activity) 
moves  the  calculated  data  closely  toward  the  experimental  curve. 


B .  High  Dose  Cases 

Figure  20  shows  the  comparison  of  calculated  and  experimental  re- 

*  10  2 

suits  for  the  case  of  the  10  ions/cm  dose  boron  implant  annealed  at 
000°C.  The  slower  electrical  activation  is  a  consequence  of  the  small 
dissolution  rate  of  boron  precipitates  and  diffusion  effects,  as  we  will 
see.  There  are  plateaus  on  the  electrically  active  boron  profiles  cal¬ 
culated  at  1,  3,  10,  and  35  minutes.  In  particular,  we  compare  the 
calculated  electrically  active  boron  profile  at  35  minutes  with  the 
corresponding  experimental  electrical  carrier  concentration  profile 
represented  by  the  squares.  The  theoretical  results  are  in  good  agree¬ 
ment  with  the  experimental  data.  1+  is  apparent  that  near  the  profile 
peak  there  is  an  immobile  and  electrically  inactive  fraction  of  boron. 
This  fraction  is  composed  of  boron  precipitates  because  the  boron  con¬ 
centration  is  in  excess  of  the  solid  solubility  limit  at  the  annealing 
temperature.  Both  theoretical  and  experimental  results  exhibit,  upon 
annealing,  a  decrease  in  the  width  of  the  immobile  fraction  near  the 
profile  peak.  This  change  of  shape  is  due  to  the  dissolution  of  boron 


In  the  1016  ions/cm^  dose  cases,  the  annealed  profiles  by  Hofker  et  al 
m  have  been  corrected  for  the  30*.  ion  yield  increased  in  the  precip¬ 
itation  area  observed  by  the  authors  . 
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B  IMPLANT  AT  900°C.  The  comparison  of  calcu¬ 
lated  and  experimental  results. 


in  this  area.  The  dissolution  takes  pla  ce  because  diffusion  has  de¬ 
pleted  the  mobile  boron,  upsetting  the  equilibrium  concentrations  of 
precipitated  and  nonprecipitated  boron.  Therefore,  the  dissolution  is 
not  only  limited  by  the  dissolution  rate  of  boron  precipitates,  but  it 
i-^  also  governed  by  the  diffusion  process.  The  relative  importance  of 
this  diffusion  can  be  assessed  from  inspection  of  Fig.  21.  In  this 
figure,  we  compare  the  total  calculated  boron  concentration  at  800  °C. 
through  lO0O°C  with  the  experimental  SIMS  profiles.  It  is  apparent 
that  at  800  °c  the  diffusion  is  very  slight  and  the  profile  peak 
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Fig.  21.  THE  ANNEALING  OF  A  10  IONS/CM  DOSE 
B  IMPLANT  IN  THE  TEMPERATURE  RANGE  OF  800 °C 
TO  1000 °C. 

remains  identical  to  the  as-implanted  profile.  At  900 °C,  the  diffusion 
is  appreciable  and  we  can  observe  a  decrease  in  the  width  of  the  top 
portion  of  the  profile  peak  upon  annealing.  At  1000°C,  the  diffusion 
is  important  and  the  removal  of  boron  from  the  high  concentration  area 
is  quite  conspicuous.  In  summary,  the  electrical  activation  of  boron 
in  the  annealing  of  high  dose  implants  is  controlled  by  both  the  dis¬ 
solution  rate  of  the  precipitates  and  by  the  diffusion  of  the  mobile 
boron . 
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C .  Low  Annealing  Temperature  Cases 

Figure  22  shows  the  comparison  of  the  calculated  and  experimental 
14  2 

results  of  the  10  ions/cm  dose  boron  implant  annealed  at  800 °C .  We 
show  the  evolution  of  the  calculated  electrically  active  boron  profiles 
at  1,  3,  10,  and  35  minutes,  and  we  compare  the  35  minute  profile  with 
the  experimental  electrical  carrier  concentration  profile  at  the  corre¬ 
sponding  time.  We  also  show  the  fit  between  the  calculated  total  boron 
concentration  profile  with  the  experimental  SIMS  profile  by  Hofker  et 
al  [2"1  .  The  agreement  is  quite  good  in  both  comparisons.  We  can  also 


14  2 

Fig.  22.  THE  ANNEALING  OF  A  10  IONS/CM  DOSE 
B  IMPLANT  AT  800  °C .  The  comparison  of  calcu¬ 
lated  and  experimental  results. 
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compare  the  electrical  activation  predicted  by  this  calculation  with  the 
experimental  activation  curve  by  Seidel  and  MacRae  in  Fig.  23.  We  can 
see  that  again  the  zero  initial  electrical  activity  in  this  calculation 
is  responsible  for  the  discrepancy  between  the  two  curves  at  short  an¬ 
nealing  times.  As  in  the  previous  instance,  we  attribute  to  room-tem¬ 
perature  annealing  the  electrical  activation  of  boron  until  the  desired 
initial  electrical  activity  is  attained,  and  we  compare  the  experimental 
annealing  curve  with  the  corrected  calculated  electrical  activation 
curve.  The  corrected  data  points  are  represented  by  'x's  and  they  are 

obtained  by  shifting  the  time  axis  by  t  =15  minutes.  We  can  see  now 

o 

that  the  corrected  calculated  points  are  in  agreement  with  the  experi¬ 
mental  curve. 
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Fig.  23.  THE  COMPARISON  OF  CALCULATED  ELECTRICAL 
ACTIVITY  VERSUS  TIME  WITH  THE  EXPERIMENTAL  ISO¬ 
THERMAL  RESULT  BY  SEIDEL  AND  MacRAE.  Dose:  1014 

ions/cm^,  800 °C . 
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Next,  we  show  in  Fig.  24  the  result  of  the  calculation  performed 

15  2 

with  the  same  set  of  parameters  for  the  10  ions/cm  dose  annealed  at 
the  same  temperature .  In  this  figure,  we  show  the  comparison  of  the 
calculated  total  boron  concentration  with  the  experimental  SIMS  profile. 
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Fig.  24.  T’E  ANNEALING  OF  A  1015  IONS/CM2  DOSE 
B  IMPLANT  IN  THE  TEMPERATURE  RANGE  OF  800 °C 
TO  1000 °C. 


We  recall  that  with  the  higher  dose  implant  the  defect  density  is  also 
higher  and  the  fraction  of  immobile  boron  trapped  by  dislocation  dipoles 
also  increases.  This  larger  fraction  of  immobile  boron  produces  a  more 
pronounced  'shoulder-like  structure'  in  the  annealed  profile.  We  also 
include  in  this  figure  the  calculated  annealed  boron  profiles  for  900°C 
and  1000 °C  to  show  that  there  is  a  significant  change  in  the  shape  of 
the  annealed  profile  as  the  annealing  temperature  is  varied.  The  agree¬ 
ment  of  the  calculated  with  the  experimental  results  is  quite  good  in 
all  cases  . 
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In  summary,  we  have  used  the  mathematical  annealing  model  in  con¬ 
junction  with  the  set  of  parameters  in  Tables  6  and  7  to  calculate  the 
annealing  of  ion-implanted  boron  under  conditions  summarized  in  Table 
3.  These  conditions  include  a  variation  of  two  orders  of  magnitude  in 
implantation  dose  and  annealing  temperature  in  the  range  of  800°C  to 
1000 °c .  In  the  calculations,  we  use  initial  conditions  and  parameters 
that  are  in  agreement  with  experiments  pertinent  to  the  annealing  prob¬ 
lem.  The  main  group  of  parameters  are  obtained  using  procedures  that 
are  based  on  thermodynamic  and  physical  considerations.  In  other  in¬ 
stances,  we  use  a  single  parameter  to  represent  the  composite  effects 
of  several  phenomena.  Naturally,  more  detailed  modelling  is  possible 
in  these  cases  but,  based  on  the  good  agreement  between  the  calculated 
and  the  experimental  results,  it  is  apparent  that  the  selected  set  of 
parameters  is  representing  indeed  the  major  elements  in  the  annealing 
problem  under  the  conditions  that  we  described. 
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Chapter  IV 


SUMMARY  AND  CONCLUSIONS 

A.  The  Three  Stream  Diffusion  Model 

In  the  first  chapter,  we  studied  the  annealing  of  ion-implanted 
boron  into  silicon  and  we  reviewed  the  relevant  experiments  that  are 
of  central  importance  for  the  annealing  problem.  With  this  background, 
we  then  discuss  a  general  diffusion  model  in  which  boron  atoms  interact 
with  positively  charged  silicon  vacancies  to  form  BV-pairs.  In  this 
interpretation ,  the  diffusion  of  boron  becomes  a  weighted  diffusion  with 
contributions  from  a  slow,  electi'ically  active  fraction  (substitutional 
boron)  and  a  fast,  electrically  inactive  fraction  (the  boron-vacancy 
pair).  Tiie  presence  of  positive  vacancies  in  the  model  has  two  impoi — 
tart  consequences.  First,  a  large  concentration  of  vacancies  in  excess 
of  the  equilibrium  concentration  level  will  modify  (by  virtue  of  the 
aforementioned  reaction)  the  relative  populations  of  the  two  forms  of 
boron,  hence  altering  the  overall  diffusion.  Second,  the  concentration 
of  positive  vacancies  in  the  silicon  lattice  is  a  function  of  the  Fermi 
level  in  the  bandgap.  This  Fermi  level  dependence  is  responsible  for 
high  concentration  anomalies  observed  in  diffusion  experiments. 

This  basic  three  stream  diffusion  model  is  capable  of  predicting 
ordinary  diffusion  of  boron,  proton  enhanced  diffusion,  and  a  number  of 
important  cases  of  annealing.  In  the  ordinary  diffusion  problem,  the 
prevalent  fraction  of  substitutional  boron  under  equilibrium  conditions 
results  in  high  electrical  activity  and  slow  diffusion.  In  the  proton 
enhanced  diffusion  problem,  the  steady  state  proton  beam  produces  an 
excess  vacancy  concentration,  which  in  turn  creates  a  large  concentra¬ 
tion  of  BV-pairs.  This  predominance  of  BV-pairs  results  then  in  low- 
electrical  activity  and  fast  diffusion.  A  somewhat  similar  situation 
is  found  in  the  typical  annealing  cases  of  ion  implanted  boron.  In 
these  cases,  the  majority  of  the  implanted  boron  ions  '  aact  with  nearby 
vacancies  and  form  BV-pairs.  Annealing  proceeds  then  with  the  conver¬ 
sion  of  a  large  population  of  BV-pairs  into  substitutional  boron.  Con¬ 
sequently  ,  during  annea 1  ing  ,  the  red istr  ibut  ion  of  boron  is  characterized 
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by  an  initial  enhancement  in  diffusion  that  reduces  to  ordinary  diffusion 
rates  as  the  electrical  activity  increases  to  near  100"'. 


B .  The  Relation  between  Diffusion  and  Electrical  Activity 

The  two  forms  of  boron  in  this  three  stream  diffusion  model  have 
dichotomous  attributes.  Namely,  substitutional  boron  is  electrically 
active  and  slow  diffusing;  on  the  other  hand,  the  boron-vacancy  pair  is 
electrically  inactive  and  fast  diffusing.  Because  of  this  dichotomy  and 
because  the  diffusion  of  the  total  boron  is  a  weighted  diffusion  of 
these  two  species,  a  relation  between  diffusion  and  electrical  activity 
ought  to  exist.  We  assume  near  100’'  electrical  activity  in  the  selection 
of  parameters  to  fit  thermal  diffusion  experiments.  This  ensures  that, 
when  the  abnormally  high  concentration  of  boron  vacancy  pairs  that  is 
produced  by  the  implantation  has  annealed,  the  enhanced  diffusion  (asso¬ 
ciated  with  the  initially  low  electrical  act.)  will  relax  to  ordinary 
diffusion  which  is  simply  the  near-thermal-equilibrium  solution  of  the 
set  of  equations.  The  analytical  analogue  of  this  discussion  is  carried 
out  in  Appendix  A.  The  boron  and  BV-pair  equations  are  combined  and  the 
predominance  of  over  is  used  to  arrive  at  a  single  approxi¬ 

mated  diffusion  equation  for  total  boron.  In  this  equation,  we  can 
identify  the  overall  diffusion  coefficient  D  with  a  product  of  D„,r 

and  the  fractional  concentration  of  BV-pairs.  Furthermore,  the  frac¬ 
tional  concentration  can  be  expressed  in  terms  of  electrical  activity. 
Hence,  we  arrive  at  a  diffusion  coefficient  (the  only  one  accessible  in 
an  actual  experiment,  hence  the  subindex),  which  is  a  function  of  the 
electrical  activity. 


C 

I)  =  1)  - —  -  =  D  (1  -  electrical  activity)  (6) 

exp  BV  C;v  +  C  BV 


For  instance,  at  the  outset  of  the  90()°C  anneal,  the  assumed  electrical 

activity  is  zero,  and  the  value  of  D  is  then  initially  equal  to  D„,,. 

exp  BV 

During  anneal,  as  the  electrical  activity  increases,  D  will  decrease. 

J  exp 

Finally,  when  the  electrical  activity  approaches  98f  D  will  ap- 

-15  2  CXP 

proach  the  value  of  1.1  10  cm  /sec,  the  ordinary  diffusion  coeffi¬ 
cient  of  boron  at  900 °C  measured  in  thermal  diffusion  experiments. 
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C .  High  Dose  and  Low  Annealing  Temperature  Anomalies 

In  the  typical  annealing  cases,  the  implantation  dose  is  moderate 
and  the  annealing  temperature  is  relatively  high.  These  conditions  en¬ 
sure  that  the  redistribution  and  electrical  activation  of  boron  evolves 
following  the  description  in  the  preceding  paragraph.  However,  when  the 
implantation  dose  is  increased  sufficiently,  the  resultant  impurity  dis¬ 
tribution  in  the  substrate  will  exceed  the  solid  solubility  limit  and  a 
precipitated  phase  will  form.  Naturally,  the  presence  of  this  immobile 
precipitated  phase  will  strongly  alter  the  redistribution  and  the  elec¬ 
trical  activation  behavior  of  boron.  A  different  situation  arises  when 
the  annealing  is  performed  at  low  temperature.  In  this  case,  many  kinds 
of  defects  form,  evolve,  and  interact  during  the  anneal.  The  annealing 
temperature  is  an  important  parameter  because  processes  have  different 
activation  energies  and  the  selection  of  the  temperature  determines  the 
duration  and  time  sequence  in  the  evolution  of  the  defects  and  the  like¬ 
lihood  of  interactions  between  them.  In  particular,  we  are  interested 
in  defects  that  can  interact  with  the  implanted  impurity.  At  800 °C,  we 
have  indications  that  the  dislocation  dipoles  that  form  during  annealing 
are  capable  of  trapping  boron  atoms  in  their  associated  strain  fields; 
hence,  rendering  the  boron  immobile  and  electrically  inactive.  For  this 
reason,  the  redistribution  profiles  have  special  characteristics  when 
the  ion  implanted  boron  is  annealed  at  low  temperature.  The  electrical 
activation  in  these  cases  is  controlled  by  the  release  of  boron  from  the 
strain  fields  as  the  dislocation  dipoles  anneal  out.  We  wish  to  extend 
the  three  stream  diffusion  model  to  cover  these  high  dose  and  low  anneal¬ 
ing  temperature  anomalies,  since  this  constitutes  a  further  test  for  the 
basic  model. 

D .  The  Model  and  the  Parameters 

The  mathematical  model  for  the  annealing  problem  is  basically  a  dif¬ 
fusion  model  extended  to  include  the  aforementioned  anomalies.  The  equa¬ 
tions  express  in  mathematical  terms  the  diffusion  and  possible  interac¬ 
tions  among  the  species.  Consequently,  the  essence  of  the  overall  model 
lies  in  the  modelling  of  each  diffusion,  each  interaction,  and  the  param¬ 
eters  associated  with  each  phenomenon.  In  this  respect,  these  elements 
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are  the  links  to  the  physical  world  and,  in  particular,  the  plausibility 
of  the  values  of  the  parameters  is  indeed  an  indicator  of  the  validity 
of  the  model.  For  instance,  the  values  of  temperature  dependent  param¬ 
eters  should  obey  simple  activation  energy  relations  and  they  should  be 
invariant  when  we  calculate  the  annealing  behavior  of  different  dose 
boron  implants  annealed  at  the  same  temperature. 

Figure  25  shows  the  temperature  dependence  of  the  parameters  in  the 
basic  diffusion  model.  The  horizontal  axis  is  linear  in  units  of  inverse 
temperature  in  degrees  Kelvin,  and  the  vertical  axis  is  a  logarithmic 
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Fig.  25.  THE  TEMPER.  VIVRE  DEPENDENCE  OF 
PARAMETERS . 


scale.  The  horizontal  arrows  adjacent  to  the  lines  give  for  each  param¬ 
eter  the  appropriate  exponent  to  be  read  on  the  logarithmic  scale.  The 
straight  lines  in  this  semi-logarithmic  representation  indicate  that  the 
temperature  dependence  of  the  parameters  follows  indeed  a  simple  activa¬ 
tion  energy  relation.  This  result  is  extremely  important,  since  it  en¬ 
ables  the  determination  of  the  set  of  parameters  once  the  temperature  is 
specified . 

The  remaining  parameters  are  given  in  Table  7.  In  some  cases,  the 
modelling  was  not  carried  to  the  extent  required  to  make  the  temperature 
dependence  apparent.  This  is  the  case  of  the  lifetime  of  vacancies,  the 
vacancy  generation  time  constant,  and  the  time  constants  associated  with 
the  precipitation  and  dissolution  of  boron.  For  instance,  the  require¬ 
ments  on  the  latter  parameters  are  for  them  to  assume  very  small  and  very 
large  values,  respectively.  In  other  cases,  the  parameters  are  chosen 
to  represent  the  composite  effects  of  very  complex  phenomena.  Such  is 
the  case  for  the  threshold  concentration  level  for  the  trapping  of  boron 
by  dislocation  dipoles  and  the  release  of  boron  from  these  defects.  In 
these  cases  ,  the  representation  is  oversimplified  and  a  more  detailed 
model  could  be  sought  if  necessary.  However,  based  on  the  quality  of 
the  calculations  in  the  preceding  chapter,  it  is  apparent  that  the  pres¬ 
ent  scheme  captures  the  key  features  In  the  low  annealing  temperature 
and  high  dose  cases  as  well. 

K.  The  Initial  Conditions 

We  wish  to  emphasize  that  the  equations  in  conjunction  with  the  pa¬ 
rameter  set  constitute  a  diffusion  model  with  provisions  for  high  dose 
and  low  annealing  temperature  anomalies.  But,  to  actually  solve  a  par¬ 
ticular  problem,  an  ordinary  diffusion,  a  proton-enhanced  diffusion,  or 
an  anneal,  we  need,  in  addition,  the  set  of  appropriate  initial  condi¬ 
tions  for  the  specific  problem.  In  the  preceding  chapter,  we  solved 
nine  cases  of  annealing  with  initial  conditions  for  boron  implanted  into 
silicon  at  room  temperature.  The  implantation  condition  is  of  impor¬ 
tance  because  the  damage  produced  by  a  low  mass  ion  is  light  and  there 
is  room  temperature  annealing  as  the  damage  is  being  produced  .  Under 
these  conditions,  it  is  possible  to  assume  that  the  initial  vacancy 
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concontra t  ion  is  in  thermal  equilibrium  in  the  silicon  substrate.  This 
assumption  will  certainly  be  inadequate  for  a  liquid  nitrogen  implant  or 
for  the  ease  of  a  heavy  ion  or  a  high  dose  implant  that  produces  an  amor¬ 
phous  damage  layer. 

F.  Discussion 

Let  us  examine  a  number  of  models  reported  in  the  literature.  These 
models  predict  different  aspects  of  the  annealing  problem  that  are  in¬ 
cluded  in  our  formulation.  Table  9  is  a  summary  of  models  proposed  for 
the  annealing  or  diffusion  of  boron  in  silicon.  We  list  the  number  of 
species  in  the  first  column.  In  the  second  and  third  columns,  we  indi¬ 
cate  whether  the  model  is  capable  of  predicting  the  redistribution  or  the 
electrical  activation  of  boron  during  the  anneal.  In  the  fourth  column, 
we  list  a  representative  source.  And,  the  last  column  contains  remarks 
pertinent  to  the  particular  model.  It  is  apparent  from  the  inspection 
oi  "ab 1 e  9,  that  the  proposed  annealing  model  predicts  both  the  redis¬ 
tribution  and  the  electrical  activation,  as  a  realistic  model  should. 

It  is  also  apparent  that  the  proposed  model  is  more  complex.  In  this 
regard,  we  wish  to  offer  ttie  following  discussion. 

Our  annealing  model  is  devoted  to  the  diffusion  of  boron  and  to  the 
annealing  of  ion  implanted  boron  under  a  wide  range  of  implantation  dose 
and  annealing  temperatures.  For  this  reason,  a  substantial  portion  of 
this  model  is  dedicated  to  the  anomalies  that  arise  when  the  implanta¬ 
tion  dose  and  annealing  temperature  vary  over  a  wide  range.  In  this  re¬ 
el'd  .  on  one  hand,  it  may  be  desirable  to  extend  the  model  to  cover  other 
anomalies .  And,  the  general  nature  of  the  numerical  implementation  of 
the  solution  enables  additions  to  be  made  with  minor  modifications  of 
t ae  computer  program.  On  the  other  hand,  deletions  could  be  appropriate 
under  other  circumstances.  For  instance,  if  the  high  dose  and  low  an¬ 
nealing  temperature  cases  arc  of  little  interest,  then  the  three  stream 
diffusion  model  with  a  subset  of  parameters  will  still  predict  ordinary 
diffusion,  proton  enhanced  diffusion,  and  the  typical  annealing  cases  of 
iioroti  implants  at  room  temperature.  This  reduced  model  is  still  very 
important  and  attractive  because,  once  the  parameters  associated  with 
the  anomalies  are  deleted,  the  remaining  parameters  have  simple  activation 
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energies  and  they  are  determined  by  the  specification  of  a  single  varia¬ 
ble:  the  temperature. 

Lastly,  we  would  like  to  elaborate  on  the  boron-BV-pa ir  reaction, 
the  principal  interaction  in  the  three  stream  diffusion  model.  In  the 
annealing  problem,  the  spontaneous  conversion  of  BV-pairs  into  substitu¬ 
tional  boron  is  driven  by  the  reaction  rate,  which  is  rather  large  dur¬ 
ing  the  anneal.  From  this  viewpoint,  it  is  apparent  that  the  boron-BV- 
pa  ir  reaction  is  the  dominant  process  in  the  annealing  problem.  However, 
in  other  problems,  especially  under  thermal  equilibrium  when  the  boron- 
BV-pair  reaction  rate  is  small,  it  may  be  appropriate  to  reexamine  the 
predominance  of  the  boron-BV-pa ir  reaction  over  other  possible  effects. 
For  instance,  depending  on  the  nature  of  the  problem  and  the  features 
one  wishes  the  model  to  predict,  it  may  be  necessary  to  consider  effects 
such  as  the  influence  of  internal  electric  fields  or  diffusion  induced 
stresses  due  to  the  mismatch  of  the  impurity  atoms  in  the  lattice. 

G .  Conclus ion 

The  annealing  model  is  capable  of  predicting  boron  impurity  profiles 
under  conditions  of  ordinary  diffusion,  proton  enhanced  diffusion,  and 
radiation  enhanced  diffusion  of  the  sort  that  occurs  when  boron  is  im¬ 
planted  at  room  temperature  and  subsequently  annealed.  In  the  latter 
case,  the  model  also  predicts  the  evolution  of  the  electrically  active 
boron  as  a  function  of  space  and  time.  In  reference  to  practical  implan¬ 
tation  and  annealing  conditions,  the  model  has  been  tested  for  all  the 
cases  in  Table  2.  In  this  wide  range  of  dose  and  annealing  temperatures, 
the  most  relevant  case  in  device  fabrication  are  in  the  lower  left  cor¬ 
ner.  The  remaining  cases  are  perhaps  less  important  but  not  less  inter¬ 
esting  since  the  abnormalities  are  a  real  challenge  to  the  annealing 
model.  In  calculating  these  cases,  we  show  that  by  extending  the  basic- 
approach  it  is  possible  to  generalize  the  model  to  predict  more  complex 
situations.  It  is  also  remarked  that  th  i  —  mode  1  and  the-  associated  pa¬ 
rameter  set  are  universal  in  the  sen^ethat  the  same'  equations  and  param¬ 
eters  are  used  in  the  prediction  of  ordinary  diffusion,  proton  enhanced 
diffusion,  and  the  annealing  of  i on- implanted  boron  into  silicon  at  room 
tempera  ture . 
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Suggestions  for  Future  Work 


It  is  apparent  that  this  annealing  model  enables  the  use  of  the 
accurately  calculated  as-implanted  profiles  in  the  calculation  of  anneal¬ 
ing  profiles  for  the  case  of  boron  implants  at  room  temperature.  In  this 
manner,  the  device  designer  can  make  a  better  use  of  the  range  calcula¬ 
tions  that  are  available  for  a  large  number  of  impurity-substrate  combi¬ 
nations.  From  this  viewpoint,  it  seems  reasonable  to  suggest  that  perhaps 
similar  annealing  models  can  be  developed  for  other  commonly  used  impur- 
it ies  . 

Phosphorus  and  arsenic  are  the  most  commonly  used  n-type  dopants  in 
silicon.  It  is  generally  accepted  that  their  diffusion  in  silicon  is 
also  governed  by  a  vacancy  mechanism.  Furthermore,  a  large  body  of  work 
on  ordinary  diffusion,  diffusion  anomalies,  lattice  location,  and  an¬ 
nealing  behavior  of  ion  implanted  phosphorus  and  arsenic  is  available  in 
the  technical  literature.  It  is  therefore  likely  that  a  careful  review 
of  relevant  experiments  will  yield  key  attributes  for  the  diffusion  and 
the  annealing  of  these  impurities.  In  this  regard,  perhaps  the  subject 
of  possible  interactions  between  impurities  and  charged  vacancies  should 
be  examined.  Then,  with  some  basic  assumptions,  the  simplest  model  with 
the  most  essential  attributes  should  be  formulated  and  solved.  Undoubt¬ 
edly,  many  iterations  of  the  described  sequence  will  be  necessary,  and 
the  availability  of  the  ’diffusion  solver’  from  this  work  will  be  use¬ 
ful  . 

We  wish  to  highlight  some  similarities  and  differences  between  the 
suggested  study  and  the  case  of  boron.  The  diffusion  of  phosphorus  and 
arsenic  exhibit  high  concentration  and  background  doping  effects.  Con¬ 
sequently,  it  is  likely  that  the  Fermi  level  dependence  comes  from  in¬ 
teractions  of  the  impurity  with  negatively  charged  vacanc ies  .  It  is  also 
possible  that  more  than  one  kind  of  impurity-vacancy  complexes  coexist 
in  the  substrate  and  their  different  diffusive  attributes  give  rise  to 
anomalies  in  the  diffusion  profile.  Certainly,  the  most  important  dif¬ 
ference  between  the  boron  and  the  phosphorus  or  arsenic  cases  i-  the 
damage  that  results  from  the  implantation.  The  heavy  phosphorus  or  ..r- 
senic  ions  produce  extended  damage  clusters  in  the  lattice  as  the  ions 
cone  to  rest.  In  these  cases,  the  room  temperature  annealing  rate  is 


91 


insufficient,  and  the  individual  damage  clusters  will  overlap,  produc¬ 
ing  an  amorphous  layer.  During  annealing,  the  implanted  impurity  will 
anneal  in  both  the  amorphous  and  crystalline  regions  following  differ¬ 
ent  behaviors.  Meanwhile,  the  amorphous  layer  will  also  anneal.  The 
regrowth  of  amorphous  silicon  layers  is  also  a  complicated  matter. 
Backsoattering  experiments  [-14]  indicate  that  the  regrowth  rate  is  de¬ 
pendent  on  the  crystallographic  direction,  and  it  is  modified  by  the 
doping  effects  of  impurities. 

Perhaps  the  problem  should  be  addressed  in  three  phases .  A  first 
phase  devoted  to  the  ordinary  diffusion  of  phosphorus  and  arsenic  and 
the  most  important  diffusion  anomalies.  A  second  phase  dedicated  to 
the  study  of  the  regrowth  of  amorphous  Si  layers.  The  outcome  of  this 
study  will  certainly  be  useful  for  the  understanding  of  the  annealing 
behavior  of  boron  that  is  implanted  at  liquid  nitrogen  temperature. 
And,  lastly,  in  the  third  phase,  the  results  of  the  first  and  second 
phases  can  be  combined  to  develop  an  annealing  model  capable  of  pre¬ 
dicting  impurity  profiles  in  the  presence  of  an  amorphous  layer  pro¬ 
duced  by  the  implantation. 
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Appendix  A 


THE  NUMERICAL  INTEGRATION  OF  THE  COUPLED  DIFFUSION  EOLATIONS 


In  Chapter  II,  we  formulated  the  annealing  model  which  consist?  of 
a  set  of  nonlinearly  coupled  partial  differential  equations.  We  will 
describe  now  the  numerical  integration  of  this  mathematical  problem.  The 
solution  consists  of  two  distinct  steps.  First,  we  transform  the  set  of 
Partial  Differential  Equations  (PDE)  into  a  larger  set  of  Ordinary  Dif¬ 
ferential  Equations  (ODE),  and  then  we  use  GEARR  [4 5 1 ,  a  FORTRAN  subrou¬ 
tine  package,  to  integrate  the  system  of  ordinary  differential  equations 


1.  The  Transformation  of  Partial  Differential  Equations  into  Ordinary 
Differential  Equations 

To  illustrate  this  transformation,  we  will  use  the  three  stream 
diffusion  model  described  in  F.q .  (Al).  This  simpler  case  is  chosen  for 
the  purpose  of  illustrating  the  problem  with  a  minimum  of  algebraic  com¬ 
plexities  . 
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In  this  system  of  PDE,  the  dependent  variables  are:  the  concentra¬ 
tions  ri  f  electrically  active  boron  (C  (x,t)),  positive  vacnacies 
( f  ( x ,  t ) )  ,  and  BV-pairs  (C  (x,t)).  Time  and  space  are  the  independent 

i>  n\ 

variables  in  this  problem.  The  transformation  under  consideration  is 
based  on  the  elimination  of  the  spatial  variable  by: 
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(a)  the  discretization  of  the  space  in  N  partitions 

(b)  the  definition  of  new  dependent  variables 

(c)  the  formulation  of  the  problem  in  terms  of  the  new 
variables 

For  example,  first,  the  x-axis  is  divided  into  33  partitions  of  width 
x.  Second,  we  define  the  concentrations  of  active  boron,  positive  va¬ 
cancies,  and  BV-pairs  in  each  one  of  the  33  partitions  as  new  variables 
Since  3  variables  are  associated  with  each  partition,  and  there  are  33 
partitions,  the  total  number  of  new  dependent  variables  is  99,  It  is 
convenient  to  define  a  vector  y  composed  of  the  new  dependent  varia¬ 
bles  ,  as  is  illustrated  in  Table  Al .  We  obtain  the  system  of  ODE  in 
terms  of  the  new  variables  by  rewriting  Eq .  (Al)  for  each  one  of  the  33 
partitions  and  by  approximating  the  second  partial  derivative  with  the 
central  differences  f  19 )  in  Kq .  (A3). 


(A3) 


In  this  equation,  the  superscripts  1,  c,  and  r  denote  left,  center, 
and  right,  respectively:  .  jc  is  the  width  of  the  partition.  The  central 
difference  approximation  in  Kcj .  (A3)  is  indicating  explicitly  that  the 
second  partial  derivative  of  the  variable  C  at  the  central  partition 
is  only  referenced  to  the  same  variable  C'  at  the  same  location  and  at 
adjacent  locations  (left  and  right).  Following  these  instructions  and 
using  Table  Al  to  express  the  new  variables  in  place  of  the  old  ones ,  we 
obtain  Eq .  (A3),  the  system  of  99  coupled  ordinary  d i f ferent ia  1  equations . 
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The  user  should  send  a  blank  magnetic  tape,  together  with  character, 
channel,  and  density  specifications,  and  ask  for  the  appropriate  ver¬ 
sion  of  GEARB.  For  a  complete  description  of  GEARB  and  the  user's 
instructions,  the  reader  is  referred  to  Ref.  45. 
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problems  in  which  the  stiffness  is  apparent,  the  use  of  the  backward 
differentiation  formula  with  the  user  supplied  Jacobian  is  more  effi¬ 
cient.  The  annealing  problem  falls  into  this  category.  In  other  prob¬ 
lems,  where  the  stiffness  is  only  apparent  in  the  later  stage  of  the 
integration,  it  is  possible  to  change  the  method  of  integration  during 
the  computation.  The  user  need  only  to  specify  the  initial  step  size 
of  the  integration.  Thereafter,  the  package  adjusts  the  step-size  ac¬ 
cording  to  the  error  incurred  in  the  integration.  The  communication 
with  GEARB  is  performed  through  the  subroutine  DRIVEB  by  means  of  a 
CALL  statement.  The  user  supplies  a  MAIN  PROGRAM  and  two  subroutines: 
DIFFUN  and  PDB .  The  main  program  controls  the  input-output,  initial¬ 
izes  variables,  selects  the  options,  and  performs  the  communication 
with  DRIVEB.  DIFFUN  defines  the  system  of  ordinary  differential  equa¬ 
tions  [Eq.  (A3)1  and  PDB  is  the  Jacobian  fEq.  (A6)l  stored  by  diagonal 
lines.  In  the  event  that  the  Jacobian  cannot  be  coded,  it  is  possible 
to  select  an  option  that  generates  the  Jacobian  internally.  The  pen¬ 
alty  associated  with  this  option  is,  of  course,  a  loss  in  efficiency. 

In  our  case,  the  sample  codes  that  we  list  in  the  next  section  do  con¬ 
tain  the  user  supplied  Jacobian. 

It  is  apparent  from  this  outline  that  the  interface  with  GEARB  has 
a  modular  structure.  The  definition  of  the  mathematical  problem  is 
coded  in  DIFFUN  and  PDB,  and  the  adaptation  of  the  mathematical  formu¬ 
lation  to  the  particular  problem  is  accomplished  in  the  main  program. 
For  these  reasons,  the  potential  user  will  find  that  DIFFUN  and  PDB  are 
of  immediate  use;  however,  lie  may  wish  to  modify  the  main  program,  to 
alter  the  width  of  the  partition,  the  step-s  i/e ,  the  error  control ,  etc  .  , 
to  suit  his  particular  problem.  In  Section  3,  we  include  the  listing 
of  the  three  stream  diffusion  model  and  the  complete  annealing  model. 
GEARB  is  available  in  both  single  and  double  precision  versions  from: 

Argonne  Code  Center 

Argonne  National  laboratory 

97ni)  South  Cass  Avenue 

Argonne,  Illinois  (>ol,39 

Telephone:  (BIB)  739-7711,  Ext.  JBGG 
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Schematically,  the  nonzero  elements  of  the  Jacobian  are  represented 
by  'x'  in  Table  A2 .  And,  the  banded  structure  of  this  matrix  is  ap¬ 
parent  . 

Table  A2 

NONZERO  ELEMENTS  OF  THE  JACOBIAN 


2 .  The  Solution  of  the  System  of  Ordinary  Differential  Equations 

The  numerical  integration  of  the  system  of  ordinary  differential 
equations  is  performed  on  a  digital  computer  using  GEARB  rial,  a  FORTRAN 
subroutine  package.  GEARB  is  composed  of  seven  subroutines  designed  to 
solve  systems  of  ordinary  differential  equations  in  which  the  Jacobian 
has  a  banded  structure.  This  structure  allows  the  storing  of  the  ele¬ 
ments  of  the  Jacobian  by  diagonal  lines,  hence  saving  storage  and  com¬ 
putation  time.  The  package  offers  several  options  with  regard  to  the 
method  of  integration.  For  instance,  the  Adams  method  with  functional 
iteration  is  recommended  for  nonstiff  problems.  On  the  other  hand,  for 


1 


0 


i  •  () 


Jacobian  of  the  system,  it  will  be  apparent  that  this  particular  order¬ 
ing  of  species  in  partitions  does  produce  a  set  of  ODE  with  a  banded 
Jacobian  in  which  the  width  o  the  band  is  small.  And,  such  a  system 
can  be  solved  very  efficiently  since  it  does  not  require  the  storage 
of  elements  outside  of  the  band. 

It  is  more  concise  and  convenient  to  use  the  vector  notation  in 
the  description  of  the  system  of  ordinary  differential  equations.  For 
instance,  Eq .  (A3)  can  be  expressed  as: 

y.  =  f.ty,.  •••,  y . . y..)  id,  ....  JO  (A4) 

ill  J  N 


where  N  is  99  in  our  example . 

The  Jacobian  of  the  set  of  ordinary  differential  equations  in  Eq  . 
(A4)  has  elements  J.  .  given  by: 


if  i 

Ji.i  ~ 

Hence , 

JH  =  -2K1  -  Ky (2 )  =  df1/dy1 
Ji2  =  _Ky(1)  = 


(A5) 


J13  "  K'  ■ 

JM  "  K1  *  'VJyi 

!  etc . 


And,  in  general,  the  nonzero  elements  of  J  are  for  i  =  (0  ...  32): 

J  (3  i  4.  1  ,3i  +  1)  =  -2K1  -  Ky(3i  +  2) 

J (3 i  +  1 , 3 i  +  2)  =  -Ky(3i  +  1)  (A6) 
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Appendix  B 


CALCULATION  OF  AN  EFFECTIVE  DIFFUSION  COEFFICIENT 


Equations  (Bl)  and  (B2)  are  equations  (Ala)  and  (Ale)  in  Appendix 

A  rewritten  in  terms  of  fluxes,  J  and  J 

B  BV 


^B  =  ^B  |  BV  -  koCBCV^ 
dt  dx  +  \  -r  j 


(Bl) 


dC  dj 

_ BV  _  BV 

( i  t  ~  d  x 


fCBV  -  Wv*' 


(B2) 


Addition  of  Eqs .  (Bl )  and  (B2)  yields: 


—  (C  +  C  )  =  —  (J  +  J  ) 
dt  B  BV  dx  B  BV 


(B3  ) 


Expressing  the  fluxes  in  terms  of  concentration  gradients  and  dif- 
fussion  coefficients,  we  obtain  Eq  .  (B4). 


•L  +  -J. 


BV 


I) 


B 


x 


+ 


D 


BV 


(BD 


For  this  equation,  we  can  show  that  the  last  term  dominates  as  follows: 

D  is  much  greater  than  I)  ,  at  000 °C;  their  values  are  7  10  ^  and 
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1.6  •'  10  “  cm- /sec,  respectively.  At  the  outset  of  the  anneal,  C 

C  and  1C  /  ix  •  >C  /  >x ,  hence  we  can  neglect  the  first  term.  As  the 
B  BV  B 

annealing  proceeds,  the  diffusion  of  boron  will  approach  ordinary  dif¬ 
fusion  conditions,  namely, 


0.03 


The  inequality  between  gradients  has  reversed;  however,  because  D  is 
much  greater  than  ,  the  last  term  still  overwhelms  the  first  one. 

We  can  then  approximate  the  total  flux  as  in  Eq .  (5): 


127 


B  BV  BV  JX 


We  now  let  P  be  the  fractional  concentration  of  BV  pairs,  which 
is  a  function  of  the  electrical  activity: 


P 


BV 


'BV 


+  C„ 


=  i  -  a 


(B6) 


Then  , 


BV 


=  =  (c 


BV 


CB> 


(B7) 


and  the  gradient  becomes: 


C. 


^  =  (c 

»x  BV 


’X 


-  V  + 


)(C  +  C  ) 

BV  B 

Ox 


CBS) 


We  can  now  identify  several  situations  in  which  the  first  term  in  Fig. 
(BS)  is  negligible  in  comparison  with  the  last  one.  For  instance,  under 
equilibrium  conditions,  we  can  rewrite  Eq .  (3)  in  Section  III  as: 


-2£  =  k  c  . 

CB  °  V 


Then,  substitution  of  the  above  equation  in  the  expression  for  1,  • 
yield s 


1 


1  + 


=  1  + 


BV 


k  +  Cv- 
o  \ 


■■lie  re 


v- 


=  V(T) 


exp 


In  general,  the  Fermi  level,  as  a  function  of  the  acceptor  concentra¬ 
tion,  will  be  a  function  of  distance;  however,  for aimea 1 ing tempera tues 


near  900 °c  and  the  boron  concentration  below  ~lo  atoms/cm  ,  this  de¬ 
pendence  is  very  weak,  hence  C^+  ~  ,  J'/Jx  ~  0,  and  the  first  torv 

in  Eq .  (B8)  is  negligible.  An  analogous  case  arises  when  the  annealing 
temperature  is  high,  above  1100°C.  In  this  case,  the  Fermi  level  is 
fixed  in  the  middle  of  the  bandgap. 

Under  nonequilibrium  conditions,  the  comparison  of  the  terms  in  Fq 
(B8)  can  be  performed  numerically.  For  the  typical  annealing  case  in 
Table  3,  the  results  of  the  calculation  show  that  the  approximation 
in  Eq .  (B8)  is  correct  under  nonoquilibrium  and  equilibrium  conditions, 

failing  only  under  equilibrium  conditions  when  the  boron  concentration 

”TS  '  3 

is  near  and  above  10  atoms/cm  .  In  other  words,  when  the  boron  con¬ 
centration  is  high  and  the  Fermi  level  is  a  function  of  distance,  !-/ 
ix  may  no  longer  be  small . 

For  cases  in  which  d?,/dx  may  be  neglected,  we  may  then  simplify 
Eq .  (B«)  to 


JC, 


BV 


<KC, 


Jx 


~  f- 


BV 


*V 


Jx 


(B9 ) 


Substitution  of  Eq .  (B9)  in  Eq .  (B5)  and  Eq .  (B5)  in  Eq  .  (B3)  then  yields 


d(cB  *  W 

dt 


=  (D 


BV 


?) 


rf2(CBV  +  V 


OX 


(BIO) 


where  (I)  •  : )  can  be  identified  with  the  overall  diffusion  coeffi¬ 

cient  I)  .  Hence, 
exp 


I) 

exp 


dbv(1  "  °° 


(Bll) 
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